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ABSTRACT 


A general  method  for  extracting  particle  size  dependent  information 
from  experimental  rate /formulation  data  was  developed  from  the  statistical 
methodology.  This  technique  was  employed  to  correlate  the  data  bases  of 
Miller.  Results  showed  thatiby  employing  an  interaction  parameter  of  4 that 
both  additive  and  additive  free  data  could  be  correlated  to  standard  error  of 
estimate  below  10.  5%.  The  effect  of  steady  radiant  energy  deposition  on 
steady  and  nonsteady  burning  was  explored.  Results  showed  that  if  the 
radiant  energy  deposited  in  the  reactive  zones  is  negligible  (an  excellent 
assumption  for  low  signature  propellants)  the  effect  of  radiant  energy 
deposition  can  always  be  accounted  for  by  substituting  the  radiation  augmented 
initial  temperature-^  T^  " T0  * ^/[pou.0}  7 "for  the  initial  propellant  tempera- 
ture T0.  1 
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FINAL  REPORT 


for 

STATISTICAL  ANALYSIS  OF  STEADY  STATE  COMBUSTION 
OF  COMPOSITE  SOLID  PROPELLANTS 


INTRODUCTION 

Background 

Solid  rocket  technology  is  in  a process  of  dynamic  development  that 
is  driven  by  current  problems.  This  process  is  particularly  evident  in 
applications  of  solid  rockets  to  tactical  weapons  because  of  recent  emphasis 
on  reduced  visible  exhaust  signature.  This  emphasis  has  led  to  elimination 
of  significant  amounts  of  condensed  phase  material  from  the  products  of 
combustion  and  the  subsequent  creation  of  a number  of  problems. 

o To  maintain  specific  impulse  parity  (or  to  minimize 
degradation)  relative  to  an  "equivalent"  metallized 
formulation  the  total  solids  content  must  be  increased. 

This  can  lead  to  physical  property  and  processing 
problems. 

o Replacement  of  metal  additive  with  oxidizer  has 

altered  the  relationship  among  rate,  formulation,  and 
environmental  variables.  Since  the  bulk  of  empirical 
knowledge  of  these  relations  resides  with  metallized 
formulations,  propellant  formulation  problems  have 
arisen. 

o Elimination  of  condensed  phase  particulates  from  the 
products  of  combustion  degrades  stability  margins 
(particularly  at  higher  frequencies)  because  particle 
damping  is  "eliminated".  As  a consequence,  com- 
bustion instability  has  become  a very  significant  de- 
sign factor. 


Clearly,  it  would  be  (and  is)  more  difficult  to  design  a solid  rocket  to  speci- 
fied constraints  with  low  signature  propellant  than  with  metallized  propellant 
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The  design  of  solid  rockets  always  strives  for  a near  optimum  for 
imposed  constraints.  The  constraints  are  set  by  the  level  of  available  technology. 
It  is  important  to  note  that  there  are  two  basic  parts  to  available  technology. 

First,  there  is  that  part  concerned  with  the  performance  limits  of  available 
material;  for  example,  the  tensile  strength  of  a case  material  and  the  theoret- 
ical specific  impulse  of  a propellant.  Second,  there  is  that  part  concerned 
with  the  way  available  materials  are  arranged  into  the  entity  we  call  a solid 
rocket  motor.  Recognition  of  these  two  factors  is  important  because  available 
information^ ) suggests  that  insofar  as  propellant  energetics  are  concerned 
the  first  path  is  peaking.  Moreover,  the  achievement  of  increased  performanc  e 
along  this  path  is  complicated  by  hazards  and  cost  consideration  (1).  Conse- 
quently, as  time  passes  superior  performance  will  become  increasingly  depen- 
dent upon  the  potential  of  inert  parts  and  design  excellence. 

These  are  not  academic  topics;  cost  robs  us  of  resources  and,  relative 
to  weaponry,  performance  can  mean  our  life. 

Recent  design  experience  has  amply  demonstrated  that  combustion 
instability  is  a major  factor  in  solid  rocket  design.  This  may  seem  surprising 
in  view  of  the  technology  that  has  been  developed  to  treat  this  problem  area  to 
date.  However,  the  simple  fact  is  that  existing  technology  has  proven  to  be  too 
cumbersome  for  low  signature  systems.  This  stems  largely  from  the  fact  that 
the  direction  of  combustion  instability  technology  was  shaped  by  instability 
problems  in  strategic  missile  system s - -that  is,  in  systems  with  highly 
metallized  formulations,  nearly  neutral  grains  ' # and  small  environmental 
temperature  ranges.  In  these  systems  instability  needs  to  be  considered  at 
only  the  lower  longitudinal  mode  frequencies  (particulate  damping  suppresses 
the  higher  mode  frequencies),  and  at  one  pressure  and  one  initial  temperature. 
Therefore,  the  number  ol  coni  binations  for  which  data  are  required  is  small. 
Taking  pressure,  initial  temperature,  and  frequency  as  variables  the  number 
of  combinations  is  on  the  order  ol  1-3.  On  the  other  hand,  in  a low  signature 
tactical  system  instability  is  not  limited  to  the  lower  frequency  modes.  More- 
over, non-neutral  traces  are  common  (boos t/ susta in),  and  the  envi ronmental 
temperature  range  is  substantial  (-/0  to  +1/0°F).  Consequently,  the  number 
of  combinations  is  on  the  order  of  20-30.  In  short,  the  designer  of  a low 
signature  tactical  system  is  looking  at  a task  that  is  roughly  an  order  of  mag- 
nitude more  involved  than  the  designer  of  a strategic  system  in  order  to  assure 
the  same  "surprise  free"  design.  Consequently,  design  procedures  that  "work 
like  a champ"  with  strategic  systems  can  simply  be  overwhelmed  by  the  compu- 
tational and  data  demands  of  a low  signature,  tac  tical  system. 


^Reference  2 presents  a case  study  of  a recent  reduced  smoke  n otor 
development  effort. 

neutral  grain  is  one  that  produces  a nominally  level  thrust  tii 
history  during  the  motor's  action  time. 
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The  point  here  is  simply  this.  As  solid  rocket  technology  shitts  to 


follow  the  dictates  of  field  experience  and  mission  analyses,  the  demands 
made  of  the  technology  also  shifts.  Sometimes  these  shifts  require  no  new 
technology  while  at  other  times  they  do.  The  present  is  one  of  the  times 
that  new  technology  is  important  because: 

o good  design  is  imperative  tor  near-optimum 
performance, 

o instability  is  a major  design  problem  in  low 
signature  systems,  and 

o existing  design  and  data -gathering  tools  are  too 
cumbersome  ior  the  iunding  levels  and  develop- 
ment schedules  ol  tactical  systems. 

To  clarity  the  latter  point,  consider  the  design  process  as  it  is 
currently  practiced.  First,  non-detailed  trade  studies  are  made  to  establish 
the  general  geometric  and  propellant  properties  required  to  meet  the  design 
constraints.  Second,  a sequence  of  detailed  trade  studies  are  made  about  one 
or  more  baseline  designs  to  establish  the  "optimal"  design  in  that  baseline 
family.  Third,  a design  is  selected  for  prototype  development.  The  detailed 
studies  include  (or  should  include)  detailed  performance,  structural  integrity, 
and  linear  stability  computations.  To  carry  out  performance  predictions  one 
needs  to  know  how  burning  rate  varies  with  pressure,  initial  temperature, 
and  crossflow.  To  perform  linear  stability  calculations  (with  existing  codes! 
one  needs  to  know  the  cavity  geometry,  the  local  mass  ettlux  from  the  burning 
surface  and  the  pressure  and  velocity  coupled  response  functions.  The  response 
functions  depend  upon  pressure,  initial  temperature,  crossflow,  and  frequency. 
Consequently,  carrying  out  the  detailed  trade  studies  requires  a substantial 
amount  of  propellant  ballistic  data,  particularly  when  it  is  recognized  that 
twenty  or  more  formulations  may  be  involved  in  the  trade  studies. 

Assume,  for  example,  that  propellant  ballistic  data  are  to  be  obtained 
at  three  pressures,  initial  temperatures,  crossflows,  and  frequent  ie^.  With 
three  replications  for  statistical  significance  there  are  81  data  bits  involved 
in  defining  burning  rate.  With  the  variable  area  T-burner  technique'  1 
employing  data  at  three  area  ratios  there  are  72cf  data  bits  for  pressure* 
coupled  response  and  an  additional  243  data  bits  for  velocity  coupled  response*  . 
This  example  serves  rather  graphically  to  illustrate  the  magnitude  of  the  non- 
steady state  ballistic  data  problem  relative  to  that  of  the  steady- state*.  When 
one  considers  further  that  thirty  tests/day  is  a very  good  rate  tor  T-burner 
testing,  that  the  cost  of  a T-burner  test  is  on  the  order  of  $50,  and  that 
twenty  or  more  formulations  are  usually  scrutinized  in  a motor  development 
program,  it  is  easy  to  see  why  stability  analyses  are  always  based  on  incom- 
plete data.  The  cost  required  for  complete  ballistic  data  is  too  large  a 
fraction  of  the  total  program  cost! 


This  example  assumes  that  flow  turning  and  veloc  ity  coupling  can  be 
unconfounded. 
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Contrast  the  above  "tactical”  example  with  a "strategic”  example 
where  one  pressure,  initial  temperature  and  three  frequencies,  crossflows 
are  involved.  Then  with  three  replications  9 data  bits  are  required  to  define 
rate,  81  data  bits  are  required  to  define  pressure  coupled  response,  and  27 
data  bits  are  required  to  define  velocity  coupled  response.  Thus,  as  noted 
before,  the  magnitude  of  the  ballistic  characterization  problem  and  hence  the 
cost  and  time  involved  is  roughly  one  order  less  than  for  the  tactical  situation. 
In  contrast,  the  funding  level  for  the  strategic  motor  development  effort  is 
roughly  an  order  of  magnitude  greater  than  for  the  tactical.  Therefore,  the 
cost  required  for  complete  ballistic  data  is  a much  more  acceptable  fraction 
(down  roughly  two  orders  of  magnitude)  of  total  program  cost'1. 

Thus,  we  are  led  to  an  interesting  conclusion.  The  current  driving 
force  for  upgrading  solid  rocket  design  techniques  stems  not  from  the  strategic 
but  from  the  tactical!  This  is  not  particularly  surprizing.  It  is  no  challenge 
to  design  a Mercedes;  it  is  a challenge  to  design  a Ford  that  is  as  good  as  or 
superior  to  a Mercedes. 

The  above  shows  rather  clearly  that  a major  weakness  in  current 
solid  rocket  design  technology  is  adequate  characterization  ot  ballistic  proper- 
ties. It  is  important  to  note  that  design  studies  are  basically  quantitative 
trading  operations.  Therefore,  if  inaccurate  ballistic  properties  are  employed, 
the  trading  operation  degenerates  to  the  qualitative  level.  This  is  adequate  for 
academic  exercises;  it  is  inadequate  lor  the  design  of  propulsion  systems  that 
will  fly*'1'.  This  problem  can  be  overcome  in  two  rather  different  ways. 

o Develop  more  effective  methods  for  defining  linear 
stability  properties  experimentally. 

o Develop  more  effective  ways  to  extrapolate  from  a 
limited  data  base. 

In  actuality  both  paths  must  be  pursued  because  copious  amounts  of  high 
quality  data  are  required  to  test  the  validity  ot  the  extrapolation  procedures. 


-The  rotating  valve  burner  under  development  at  CSD(4)  offers  roughly  a 
three  fold  reduction  in  the  cost  of  characterizing  nonsteady  ballistic 
properties. 

-'-The  knowledgeable  reader  will  certainly  note  that  we  have  and  are  getting 
the  design  job  done  without  new  design  tools.  However,  it  is  not  without 
considerable  "cut  and  try"  at  the  prototype  motor  level.  This  is  expensive. 
Moreover,  it  shifts  emphasis  away  from  an  optimal  solution  and  toward 
any  solution.  This  writer  believes  that  the  most  critical  aspect  of  the 
design  process  is  the  detailed  trade  studies  and  that  realism  in  this  phase 
is  critical  to  the  outcome. 


AFOSR  Statistical  Combustion  Modeling  Program 


This  program  has  been  aimed  at  steady-state  combustion  of  poly- 
disperse  heterogeneous  propellants.  The  general  methodology  has  been  to 
employ  a statistical  approach  to  relate  the  areal  mean  burning  rate  of  the 
propellant  to  the  areal  mean  burning  rate  of  monodisperse  pseudo-propellants 
whose  properties  are  derived  from  the  statistical  formulation.  The  significance 
of  the  arrangement  is  that  it  is  much  easier  to  model  propellant  with  a single 
geometric  parameter  than  propellant  with  a distribution  of  characte ristic 
dimensions.  A monodisperse  propellant  combustion  model  has  been  constructed 
from  the  BDP  mod  el  and  combined  with  the  aior  ementioned  statistical  pro- 
cedure to  yield  a steady-state  combustion  model  for  polydisper se,  AP-hydro- 
carbon  binder  composite  propellants.  The  model  has  been  tested  against  the 
extensive  data  bases  generated  by  Miller,  et.  al.  ^).  it  has  been  found  that 
the  model  is  capable  of  quantitative  predictions  of  the  effect  of  formulation 
variables  on  both  rate  and  exponent  for  additive  free  formulations.  However, 
predictions /cor relations  of  the  data  bases  with  additives  (aluminum,  iron  oxide} 
were  generally  poor^)#  In  addition,  it  should  be  noted  that  formulations  with 
a large(coarse  diameter ) /(fine  diameter)  ratio  generally  showed  poorer 
correlation.  These  defects  have  generally  been  attributed  to  the  fact  that  the 
various  pseudo-propellants  interact  more  strongly  than  presently  accounted 
for  and  that  additives  influence  these  interactions. 

Subsequent  theoretical  developments  have  attempted  to  extend  the 
theory  to  include  erosive  burning^  and  pressure^)  and  velocity  coupled 
nonsteady  burning  ' Insufficient  data  exists  at  present  for  detinitive  com- 
ments at  this  time.  However,  the  outlook  is  not  promising  for  the  nonsteady 
extensions.  A primary  reason  for  this  is  the  current  inability  to  come  to 
grips  with  the  nonsteady  temperature  field  in  the  condensed  phase  of  a com- 
posite propellant  without  resort  to  postulates^ ^), 

In  the  Background  section  of  the  INTRODUCTION  it  was  pointed  out 
that  a major  problem  in  the  solid  rocket  design  process  is  collection  of 
adequate  nonsteady  data.  Therefore,  the  question  ' why  pursue  steady-state 
modeling?  " arises  quite  naturally.  The  answer  is  that  potential  exists  for 
"transforming"  steady-state  data  into  nonsteady  state  "data"  for  those  situa- 
tions where  the  characteristic  time  of  the  unsteady  environment  is  large  com- 
pared to  the  characteristic  times  of  the  reactive  zones  of  the  proc  ess  . 

That  is,  a t rans  formation  should  be  possible  for  those  frequencies  where  the 


*An  areal  mean  is  defined  as  f-Um  ^ 5 fdS/Srv 

Sq— So 

- Reference  8 presents  an  excellent  summary  of  these  developments. 
::':  : See  Reference  11  for  a lucid  exposition  ol  this. 
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reactive  regions  behave  quasi- steadily  and  the  transient  aspects  are  confined 
to  the  nonreactive  condensed  phase.  This  transformation  already  exists  ior 
homogeneous  propellants^13*.  The  difficulty  is  to  extend  this  methodology  to 
a composite  propellant.  This  is,  to  all  appearances,  a formidable  problem 
because  at  any  instant  of  time  the  burning  surface  is  composed  of  all  states 
of  all  pseudo-propellants.  Therefore,  some  sort  of  multidimensional  solu- 
tion would  appear  necessary.  Unfortunately,  this  is  beyond  the  capabilities 
of  present  computational  machinery.  However,  during  the  19/6  work  a 
promising  new  approach  was  conceived.  For  steady-state  conditions,  the 
statistics  governing  an  areal  mean  at  fixed  time  are  the  same  as  those  for 
a temporal  mean1'  at  fixed  location  (ergodic  theorem(^)).  Consequently,  for 
small  deviations  from  steady- state,  as  might  occur  for  linear,  nonsteady 
process,  it  should  be  possible  to  replace  an  areal  mean  (as  employed  in  the 
steady-state  modeling)  with  a temporal  mean.  This  may  seem  of  small  con- 
sequence but  it  opens  the  door  to  a one -dimensional  methodology  for  computing 
nonsteady  response  for  composite  propellants.  Since  the  areal  mean  statistics 
are  known  for  the  pseudo-propellants,  the  probability  of  finding  any  pseudo- 
propellant in  a vertical  stack  of  pseudo-propellants  (a  'Dagwood  sandwich’  ) 
is  known.  In  addition,  for  each  pseudo-propellant  the  mean  quasi-steady 
behavior  of  the  reactive  zone  is  known  from  the  steady-state  calculations. 
Moreover,  the  Z-N  methodology( 1 ^ ) tells  how  to  carry  these  over  to  the 
nonsteady  state.  Therefore,  by  averaging  the  response  of  the  vertical  stack 
of  pseudo-propellants  to  fluctuating  external  conditions  the  temporal  mean 
is  achieved.  This  should  be  the  areal  mean  response  function  desired. 

Steady-state  combustion  modeling  is  crucial  to  this  enterprise.  To 
apply  Z-N  methodology  one  needs  to  know  the  behavior  of  the  reactive  zone 
in  detail.  Quite  frankly,  the  steady- state  statistical  combustion  model 
supplies  exactly  this  information.  Therefore,  to  the  frequency  limitation 
mentioned  previously,  steady  and  nonsteady  state  combustion  are  rather 
imtimately  connected  parts  of  the  same  phenomena. 

Hopefully,  the  reader  feels  as  enthusiastic  as  this  writer  at  this 
point.  Unfortunately  ^ there  is  a fly  in  the  ointment:  the  theoretical  combus- 
tion model  works  only  for  additive-free  formulations  at  present,  whereas 
rocket  motor  propellants  invariably  contain  additives.  Thus,  one  might 
conclude  prematurely  that  the  aforementioned  strategy  is  fit  solely  for 
academic  purposes. 


A temporal  mean  f=lim 

A*  — *oo 


t+AJ 

^ fdt /At 
t 


w i 
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As  noted  previously,  a statistical  combustion  model  consists  ot  two 
parts:  a statistical  framework  relating  pseudo-propellant  properties  to  pro- 
pellant properties  and  a combustion  model  for  computing  pseudo-propellant 
properties.  The  primary  difficulty  with  current  combustion  models  is  an 
inability  to  come  to  quantitative  grips  with  additives.  To  circumvent  this 
difficulty,  note  that  the  required  pseudo-propellant  information  must  exist 
within  an  adequate  steady-state  data  base.  Moreover,  by  treating  the  pro- 
pellant data  as  knowns  and  the  pseudo-propellant  properties  as  unknowns  the 
statistical  framework  provides  a means  for  computing  pseudo-propellant 
properties  from  an  adequate  steady-state  data  base.  This  procedure,  for 
restricted  conditions^ was  developed  in  1977  and  was  found  to  yield  quantitative 
results  for  additive  free  formulations  and  qualitative  results  for  formulations 
with  additives'^). 

Figure  1 illustrates  the  general  strategy  at  the  start  of  the  19/8 
program.  There  are  two  major  paths  to  the  goal  of  generalizing  experimental 
steady-state  data  so  that  it  can  be  employed  to  predict  steady  and  nonsteady 
properties  of  propellants  in  the  same  formulation  family  as  the  data  base. 

The  path  through  the  theoretical  model  requires  a substantially  smaller  data 
base  than  that  through  the  statistical  framework.  This  is  its  major  advantage. 

The  tinal  path  to  nonsteady  properties  is  not  operational.  The  strategy  in 
1978  was  to  push  development  along  both  paths.  However,  because  of  difficulties 
encountered  along  the  statistical  framework  path,  little  was  accom plished  on 
the  statistical  combustion  model  path. 

In  addition  to  these  tasks  work  was  undertaken  to  determine  the  effect 
of  thermal  radiation  on  nonsteady  pressure  coupled  burning.  This  is  an 
obvious  scaling  factor  that  must  be  accounted  for  in  the  application  of  response 
functions  measured  in  burners  of  small  geometric  scale  to  motors  with 
significantly  larger  geometric  scale.  In  addition  to  the  obvious  technical 
benefits  of  pinning  down  these  effects  there  is  the  attendant  benefit  of  detailed 
familiarization  with  the  Z-N  methodology  to  be  employed  in  the  subsequent 
attempt  to  complete  the  pseudo-propellant  properties  to  nonsteady  properties 
link. 

« 

Accomplishments  in  1978  are  described  in  the  following  section 
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Figure  1.  Schematic  Illustrating  Strategies  to  Generalize  Steady -State  Data 
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ACCOMPLISHMENTS 


I.  The  Effect  of  Interactions  in  Statistical  Combustion  Modeling 

The  design  of  solid  propellant  rockets  to  maxirrize  ’ performance' 
while  satisfying  envelope,  stability,  processing,  signature,  and  cost  con- 
straints is  a complex  task.  Information  central  to  the  successful  completion 
of  this  task  are  functional  relationships  between  the  steady  and  nonsteady 
ballistic  parameters  needed  to  predict  performance  and  stability  and  their 
independent  variables  (propellant  formulation,  pressure,  initial  propellant 
temperature,  crossflow,  frequency).  With  the  large  number  of  independent 
variables  involved  experimental  definition  over  the  variable  ranges  involved 
in  many  design  situations  are  prohibitively  expensive  with  current  techniques. 
Moreover,  the  very  nature  of  experimental  characterization  always  leads  to 
discrete  rather  than  continuous  information.  Since  information  in  discrete 
form  is  not  generally  compatible  with  optimization  strategies,  experimental 
data  must  first  be  transformed  into  a continuous  form.  Generally  speaking, 
if  all  other  things  are  equal,  a continuous  form  compatible  with  physical 
principles  is  to  be  preferred  to  those  devoid  of  insight  because  they  offer 
greater  potential  for  accurate  inte rpolation / extrapolation. 

In  addition  to  simply  correlating  experimental  data  into  a continuous 
form  a proper  theory  also  offers  potential  for  an  even  more  desirable  goal: 
transformation  of  steady-state  data  into  nonsteady  state  data.  This  goal  is 
desirable  because  nonsteady  ballistic  data  are  much  more  expensive  to 
acquire  than  steady-state  data. 

In  reference  7 a technique  was  developed  that  could  extract  ballistic 
data  in  a fundamental  form  that  readily  accounted  for  the  effect  of  particle  size 
on  steady-state  ballistic  properties  and  held  potential  for  making  the  steady 
to  nonsteady  state  transformation  possible.  Correlations  of  Miller's  additive 
free  data  base  revealed  the  ability  to  accurately  correlate  rate  and  exponent 
data.  However,  correlation  of  the  additive  data  showed  that  the  methodology 
did  not,  in  general,  work  with  propellants  containing  additives.  The  objective 
of  this  work  is  to  approach  this  problem  again  with  a more  general  version  of 
the  basic  theory. 


11 


From  statistical  combustion  theory(^)  the  mean  burning  rate  of  a 
heterogeneous,  propellant  with  polydisperse  oxidizer  is  related  to  the  burning 
rate  of  a sequence  of  monodisperse  pseudo-propellants  by 


(1-1  ) 


where 


— * 


° the  mean  burning  rate  of  the  monodisperse 

pseudo -propellant  with  oxidizer  having  0 * D * 


a/ 

° *s  ^e  oxidizer  mass  fraction  of  the  monodisperse 

pseudo-propellant  with  oxidizer  having  , and 


o dur^  is  the  mass  fraction  of  oxidizer  with 

For  most  oxidizer  grinds  it  has  been  shown  that  a log  normal  distribution 
closely  approximates  the  real  distribution^ 1 5 ).  Therefore,  to  a good  approxima- 
tion 


2 —— ==;  exp[“L(^>u  b - /rr^]/z|  diub 


M 

N 

(14) 

It  has  also  been  shown  that 

frf\  - ^ 


K%^  ' [u 4 ^ /(po] 


(1-2) 


(1-3) 


where 


(1-4) 
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Eq.  1-4  can  be  written  as 


Since  ' ~ 3fM  * \ = p*.  ( 'O  / fr.  , 


Therefore, 


r*  - f C H +[('-*«v  ? ^ 3^1  K"a  v 

C.  b (J 

Variations  in  the  environmental  variables  do  not  change  the  propellant 
recipe.  Therefore,  appropriate  differentiations  of  Eq.  (1-6)  yield 


{^P'j  --  a; » V V 


It  '•  ^ ] \SX^ 


* ■£>  i c>  b 


(I- 7a) 


( I - 7 b > 


( I - 7 c ) 


^ ( I- 7d) 


It  should  be  noted  that  Eqs.  (1-6)  and  (I-7a,b)  are,  since  they  arise  largely 
from  conservation  of  mass,  and  are  concerned  with  means,  loosely  tied  to 
propellant  structure.  On  the  other  hand,  Eqs.  (I-7c,d)  are,  since  time  is 
involved  explicitly,  more  closely  tied  to  propellant  structure.  If  the 
structure  is  random,  Eqs.  (I-7c,d)  should  hold.  If  it  is  ordered,  they 
should  not;  layer  frequency  "resonances"  would  occur. 
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Equations  1-6  - 7d  have  been  employed,  with  theoretical  combustion 
models  to  make  "a  priori"  predictions  of  propellant  ballistic  properties^'  ^ ^ 1 

Generally,  results  show  excellent  correlations  for  and  n^.  with  additive  free 
formulations  and  poor  correlations  for  formulations  with  additives.  Theory/ 
experiment  comparisons  for  cr^  t , and  R^  t are  inhibited  by  the 

inadequate  data  base.  In  addition,  best  correlations  of  theory/experiment 
occur  when  m = 3 (see  Eqs.  I -6  and  1-7).  It  is  important  to  note  that  selection 
of  m = 3 is  not  based  on  definitive  studies;  the  computational  burden  is  punitive. 


Examination  of  Eq.  (1-6)  shows  that  when  m = 3 it  becomes  the  simple 
linear  relation 


(1-8) 


where  ^is  the  mass  fraction  of  the  ktb  mode  oxidizer  in  the  formula- 

tion and  r£  is  the  burning  rate  of  pseudo-propellant  formed  from  the  ktb 
modes  oxidizer.  With  Eq.  (1-8)  simple  linear  relations  for  n^  ^ c etc 
that  involve  the  modal  pseudo-propellant  properties  follow  easily^. 
Correlations  of  experimental  rate  and  exponent  data  have  shown  that  this 
approach  worked  extremely  well  for  additive  free  formulations  and  some 
formulations  with  additives^).  The  question  to  be  asked  is  "would  m£3 
permit  better  correlation?  " If  one  attempts  to  answer  this  with  the 
theoretical  approach,  results  are  confounded  with  the  combustion  model's 
infidelities. 


Theoretical  computations ^ ^ ^ and  experimental  reductions  ^ ^ both 
suggest  that  rp  ^ is  a smooth  function  of  In  D.  Therefore,  it  is  expected 
that  a "low"  order  approximation  of  r^  ^ = f(lnD)  over  the  interval 

- 3cr  sixjbsSnub  t W would  be  adequate  for  the  task  at  hand.  The  first 

1 1 M M 1 

approach  employed  a power  series  approximation. 


M 

- >-  0.^1)  (I_< 

for  ?p # cj.  The  OL^  were  sought  with  a nonlinear  optimization  scheme 
(PATSH).  A major  difficulty  with  this  approach  became  apparent  after 
coding  was  completed;  "optimal"  correlations  gave  r'p  ^ < O for  some  D. 
Since  this  is  physically  impossible,  means  for  "forbidding"  these  "solutions' 
are  necessary.  We  could  see  no  simple1'  way  to  accomplish  this.  Conse- 
quently, this  approach  was  abandoned  in  favor  of  a piecewise  linear  approach. 
In  this  approach  the  computation  process  defines  ^*p  cj  - r\  at  specified 
In  D-  for  i = 1,  M . Therefore,  extraction  of  optimal  rj‘  can  be  easily 

guaranteed  by  employing  - \r*\  in  the  computational  process.  This  "folds 
the  r\  <o  domain  back  over  the  q >0  domain  thereby  enabling  PATSH  to 
search  over  both  positive  and  negative  r.  without  an  r*  0 constraint. 
Experience  has  shown  that  PATSH  can  become  very  troublesome  when  faced 


-Easy  to  implement  and  computationally  rapid  and  effective. 
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with  constraints  of  this  type.  Moreover,  the  extremely  valuable  feature 
of  the  polynomial  approach- -the  ability  to  evaluate  integrals  involving  D once 
and  for  all  for  each  m--is  retained.  In  addition,  once  an  approximation  for 


rp  d = f(lnD)  has  been  extracted,  the  xif  i = 1,  M can  be  relocated--to 


provide  a better  approximation  to  those  regions  where  ^ changes  rapidly 
with  In  D. 


With  the  piecewise  linear  approach  (x  = In  D) 


V r* + U*. 


(1-10) 


Therefore,  Eq.  1-6  becomes 


cAur  t 

Q 


v B*- * fel  W*  r 

b ^ u 

*•  r*  r*  1 

. rL+i  ■ t l i 

‘HI 


(1-11) 


1=1 


With  dwd  given  by  Eq.  1-2  this  becomes 


K 


M-i 


N\ 


*»■  *UH  & ^ 


(1-12) 


J x - 


Examination  of  Eq.  1-12  shows  that  there  are  four  types  of  integrals  to 
be  evaluated 


*vH 

X ■ 

(1-13) 

iL> 

*• 

(1-14) 

c 

V. 

(1-15) 

c 

^ |uT  ( * e*P  fc">  ‘ ' i d * 

(1-16) 
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Let  UL  - 


o:\rr 


then  fC  0^  4 UL  - A* 


and 


\,+\ 


^ j_  f -a 

Li.,V.  = \fir  \ e ^ 


The  error  function  is 


i f 

ert  9 yfc)*-  <H 


Thus , 


C “ £ ["^V,  ' 'rtM 


Rewrite  Eq.  1-14  as 


M \ 


liy  Vfilr  | |A‘  * \ \ e>pt'U  V^fj  ^ d-20) 


With  the  above  change  of  variable  this  becomes 


7. 

U.I 


T<*1  - r0)  \ ~i.*i  > 

Ll  ,s(e  -£  ) 


Rewrite  Eq.  1-16  as 
K’v*» 


it 

v'{5r 

* *1 

- T“' 

The  latter 

is  recognized  as 

w = expi(ft^-^)x.]/\TTf  and  dvr  = 

r *-*fc 


XUi 

* f 


Then 


i>  --  -Tfc  e 


/vr 


(I-  17) 


(1-18) 


(1-19) 


(1-21) 


(1-22) 


(1-23) 
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Consequently 


*v*» 

Vi?  SCx'^e^'3>'UH-;')Vx  '■  h e»pt('m’*)*-i('iiv)] 


Thus , 


T(4*  r-  , i,rUi  °fc.  r . / )c-tfe\V 

h,  t-  tV^-^fcJIiA  ' 4-  “*(M«  "Ur"  )] 


VM 


With  Eqs.  1-2  and  1-15 
M U) 


^H.rzvkt:-|vJx; 

b 

Attempts  to  reduce  the  remaining  integrals  to  defined  functions 
fail  because  integrals  involving  the  integral  of  the  error  function  appear. 
Consequently,  l[^  is  evaluated  by  numerical  quadrature. 

With  these  terms  Eq.  1-12  can  be  rewritten  as 

- rl/  * Tu\ 

r*  “ 2,  \(ri  -*  1 * 

r*  _ p*  tr  r - M>  at  '-“o,  . TtSM) 

+ irrr  sVl 

For  each  of  Miller’s  formulation  series  0 «©E  a-  are  fixed.  Therefore, 

can  be  evaluated  once  and  for  ail  and  ij^  can  be  computed  once  for  each 


M 


(1-25) 


(1-26) 


(1-27) 


lil) 

rn 


Equation  1-27  gives  r^  in  terms  of  the  r . . To  evaluate  the  r-  a 
nonlinear  optimizer  (PATSH)  and  Miller’s  data  base  were  employed. 
Basically,  a subset  of  Miller's  data  base  was  input,  m was  specified,  and 
the  optimizer  selected  the  r * such  that  the  standard  error  of  estimate 
between  prediction  and  data  was  minimized.  Appendix  B tabulates  the 
FORTRAN  IV  code  that  was  employed. 

Table  1-1  presents  the  standard  error  of  estimates  achieved.  The 
most  striking  feature  of  these  results  is  that  the  correlation  improves  as  m 
increases.  On  theoretical  grounds  one  would  expect  to  see  m~2;  on  the 
basis  of  limited  model  computations  it  was  found  that  m~3  effected  better 
correlation  than  m~2  (no  computations  with  m>  3 were  made);  these 
results  show  that  best  correlation  is  achieved  with  m ~4.  Note  that  with 
m = 4 that  all  of  Miller’s  rate  data  is  correlated  with  standard  error  ot 
estimates  below  10. 
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TABLE  1-1 


CORRELATION  RESULTS 


DATA  SET 

PRESSURE 

psi 

m 

SEE,  % 

COMMENTS 

SD-I-88 

500 

2.  0 

6.  34 

88%  total  solids. 

1000 

2.  0 

12.  65 

18%  24/*  Al. 

2000 

2.  0 

20.  36 

1000 

0.  5 

14.  54 

1000 

1.  0 

14.  29 

1000 

1.  5 

13.  97 

1000 

2.  5 

1 1.  38 

1000 

3.  0 

8.  13 

1000 

4.  0 

5.  48 

SD-III-88 

500 

2.  0 

4.  31 

88%  total  solids. 

1000 

2.  0 

6.  65 

no  additives 

2000 

2.  0 

1 1.  25 

500 

3.  0 

3.  05 

1000 

3.  0 

4.  49 

2000 

3.  0 

7.  57 

500 

4.  0 

2.  6 

500 

4.  5 

2.  16 

SD-IV  - 88 

500 

4.  0 

4.  03 

88"  total  solids, 

1000 

4.  0 

6.  54 

18%,  90/*  Al. 

2000 

4.  0 

10.  49 

SD-V-88 

500 

4.  0 

1.  68 

88%  total  solids, 

1000 

4.  0 

2.  36 

18%  6fi  Al. 

2000 

4.  0 

4.  52 

SD-  VI-  90 

500 

4.  0 

1.  08 

90%  total  solids, 

1000 

4.  0 

2.  53 

21"'  24 /*  Al. 

2000 

4.  0 

5.  52 

SD-  VII-  88 

500 

4.  0 

1.  39 

88"'  total  solids. 

1000 

4.  0 

3.  62 

18"  24/*  Al. 

2000 

4.  0 

3.  91 

+ 1 F e 2 O 3 

Standard  Error  of  Estimate 
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Examination  of  the  results  for  the  SD-I  and  SD-III  series  suggests 
that  correlation  is  not  strongly  dependent  upon  m for  the  additive  free 
formulations.  However,  for  formulations  with  additives  use  of  m 
substantially  improves  the  correlation  . It  is  interesting  to  note  that  in 
every  case  the  correlation  degrades  as  pressure  increases. 

Figures  1-1,  1-2,  and  1-3  present  the  computed  pseudo-propellant 
rates  as  a function  of  oxidizer  particle  size.  For  m~  3^7)  these  rates 
tended  toward  an  asymptotic  limit  as  diameter  approached  zero  and  to  zero 
as  diameter  increased.  These  results  for  m = 4 show  a much  more  com- 
plex multi- extremum  behavior.  For  the  formulations  with  aluminum  as 
the  sole  additive  burning  rate  decreases  toward  zero  with  decreasing  diam- 
eter for  diameters  below  1 /*.  However,  with  the  addition  of  1%  Fe^O^  this 
trend  is  reversed  and  burning  rate  increases  with  decreasing  diameter  for 
diameters  below  If*. 

These  calculations  have  shown  that  by  employing  m ^4  that  rate/ 
particle  size  data  can  be  correlated  and  diameter  dependent  pseudo-propellant 
properties  extracted  for  formulations  with  and  without  additives.  It  this 
technique  remains  valid  for  temperature  sensitivity  data  (insufficient  data 
currently  exists  to  say  either  yea  or  nay),  the  pseudo-propellant  properties 
needed  for  a composite  propellant  Z-N  model  can  be  extracted. 


i 
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igure  1-1.  Pseudo-Propellant  Burning  Rates  For  Miller's  SD-IV-88  Formulation  Series  When 


Pseudo-Propellant  Burning  Rates  for  Miller's  SD-V-88  Formulation  Series  When 


SD-VII-88  Formulation  Series  When 


II.  The  Effect  of  Thermal  Radiation  on  Pressure  Coupled  Response 


1 


Methodology  for  computing  the  linear  stability  margin  of  a solid 
propellant  rocket  motor  includes  testing  propellant  in  a small  scale  labora- 
tory burner  (T-burner,  SEV-burner,  rotating  valve  burner)  to  determine 
its  response  to  pressure  and  crossflow  oscillations.  These  response  func- 
tions are  then  employed  in  linear  stability  calculations  for  full  scale  motors. 
For  this  process  to  be  useful  the  response  function  in  the  motor  environment 
must  be  either  that  of  the  laboratory  burner  or  scalable  from  it.  At  present, 
little  is  known  about  this  scaling  process  because  neither  response  functions 
nor  motor  stability  margins  can  be  defined  with  precision.  Calculations  have 
demonstrated  that  the  time  mean  flow  field  in  a rocket  motor  is  dependent 
upon  its  geometric  scale*16'  and  that  crossflow  can  effect  pressure  and 
velocity  coupled  response  in  ’'strange”  ways*17'.  Therefore,  potential  exists 
tor  significant  scaling  effects.  However,  even  in  the  situation  of  pressure 
coupling  without  crossflow  there  is  still  an  environmental  difference  between 
burner  and  motor.  In  the  motor  radiosity  at  the  burning  surface  will  be 
larger  than  in  the  burner  because  the  beam  length  is  larger  and  the  "walls" 
are  hotter.  This  will  be  particularly  true  for  reduced  and  minimum  smoke 
propellants  because  their  effective  gas  emissivities  will  be  less  than  those 
of  propellants  whose  products  contain  significant  amounts  of  condensed  phase 
particulates. 

It  is  known  that  radiosity  level  can  alter  burning  rate.  In  most  cases 
this  effect  is  small  ( <10%).  Therefore,  on  these  grounds  one  might  dismiss 
the  effect.  However,  radiation,  since  it  deposits  energy  in  the  condensed 
phase,  directly  alters  the  subsurface  thermal  field.  This  may  alter  the 
frequency  dependent  character  of  the  response  functions.  Two  extreme 
situations  appear  to  exist.  If  the  extinction  depth  of  the  radiation  is  much 
larger  than  the  thermal  wave  thickness  ( the  majority  of  the  energy 

deposition  will  occur  beyond  the  thermal  wave,  the  important  cha  racte  ristic 
length  is  still  the  thermal  wave  thickness,  and  the  propellant  will  appear  to 
be  preheated.  For  a specific  "homogeneous"  propellant/1^' 


Therefore,  for  this  situation  the  effect  of  radiosity  would  he  to  alter  the  response 
function  by 

where  ST0*^  . However,  if  the  extinction  depth  is  smaller  than  the  thermal 

wave  thickness  ( &,<  X/u*  ),  the  majority  of  the  energy  deposition  will  occur 
within  the  thermal  wave;  there  will  be  two  characteristic  thermal  lengths 
i and  k/u°  );  and  possibility  for  "resonance  amplification"  exists. 
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If  propellant  is  homogeneous,  radiant  energy  will  be  deposited  in 
exponential  fashion  and  penetration  will  depend  solely  upon  the  extinction 
length  of  the  propellant.  However,  if  the  propellant  is  heterogeneous, 
radiant  energy  deposition  will  depend  upon  both  the  constituents  trans- 
missivities and  the  propellant's  structure.  The  problem  of  radiant  transfer 
in  a heterogeneous  media  is  difficult.  The  purpose  here  is  not  to  solve  this 
problem.  Rather,  the  intent  is  to  examine  qualitative  aspects  of  the  effects 
of  radiation  on  pressure  coupled  response.  Toward  this  end  two  limiting 
cases  will  be  examined.  In  the  first  the  propellant  is  considered  homogeneous 
for  both  radiant  energy  deposition  and  conductive  transport.  In  the  second, 
the  propellant  is  considered  to  be  black,  opaque  binder  and  transparent 
oxidizer  for  radiant  energy  deposition  and  homogeneous  for  conductive  trans- 
port. These  represent  simplistic  models  for  homogeneous  and  heterogeneous 
propellants. 

The  path  to  a solution  of  this  problem  will  be  to  define  radiant  energy 
deposition  in  both  cases.  For  the  heterogeneous  propellant  model  statistical 
combustion  modeling  results^^)  will  be  employed.  The  effect  of  this  radiant 
energy  deposition  on  pressure  coupled  response  wTill  be  determined  by  employing 
Z-N  methodology^  This  is  founded  upon  the  following  basic  assumptions. 

o The  propellant  is  homogeneous. 

o The  reactive  zones  behave  quasi-steadily. 

o Functions  u.  (^p;T0)  and  ^(£>*0  are  defined 
by  experimental  data. 

The  first  assumption  is  common  to  all  existing  "exact"  analyses  of  nonsteady 
combustion.  The  second  implies  both  an  upper  bound  on  frequency  for  validity 
and  that  condensed  phase  reactions  are  constrained  to  the  surface.  Relative 
to  the  latter  assumption  it  is  important  to  note  that  only  partial  derivatives 
of  these  functions  appear  explicitly.  Therefore,  these  data  must  be  very 
accurate.  Since  Tg  data  are  virtually  unobtainable  (for  heterogeneous 
propellants  the  concept  of  ji  surface  temperature  is  incorrect  and  anything 
but  accurate  and  n°  data  are  not  very  precise,  Z-N  methodology  is  currently 
inoperable  in  its  originators  context.  Consequently,  "why  pursue  the  Z-N 
approach?"  There  are  four  parts  to  the  answer.  First,  the  Z-N  method 
does  not  imply  any  reactive  zone  model.  Therefore,  it  can  employ  all. 

Second,  a path  to  predicting  the  nonsteady  response  of  heterogeneous  pro- 
pellants has  been  devised  that  employs  Z-N  concepts^  ^ ^ Third,  better 
methods  are  being  devised  to  measure  u.° T0 ) . Fourth,  the  possibility  of 
replacing  Tg  data  with  more  readily  obtainable  nonsteady  data  exists  There- 
fore, future  promise  justifies  use  of  the  Z-N  method. 


j 


14 


For  homogeneous  propellant 


d T » T A*  /CR 

Integration  and  application  of  the  boundary  condition 


~ ^ ( IF  - 2 ) 

gives 

1 ^ e ( II- 3) 


The  radiant  energy  deposition  per  unit  volume  is,  from  a radiant  energy 
balance , 


<|>  - d J ,'cU 


(II  — i ) 


With  (II- 3)  this  becomes 

$ 3 ^ & 


(11-5) 


Thus,  for  homogeneous  propellant  energy  is  deposited  in  exponential 
fashion.  Note  that  unsteadiness  in  $ is  associated  with  unsteadiness  in 
J . 

A' 

If  the  propellant  is  heterogeneous  with  black,  opaque  binder  and 
transparent  oxidizer,  radiant  energy  will  be  deposited  at  the  binder  surfac  e 
and  at  oxidizer  binder  interfaces  beneath  that  surface.  With  this  simplistic 
model  radiant  energy  deposition  will  depend  solely  upon  the  oxidizer  particle 
size  distribution. 


Consider  the  single  oxidizer  particle  illustrated  by  Figure  II- l. 
Radiant  energy  enters  the  exposed,  convex  surface  (a,  b).  It  is  assumed  that 
this  surface  is  "rough"  so  that  it  appears  to  be  a diffuse  emitter.  Thereto  re, 
the  radiant  energy  crossing  S(o)  is  also  the  energy  entering  the  subsurface. 
Since  S(o)  cannot  see  itself 


F ■*  F - l 

SC*)  ♦ ^3 


and 


S(o)  -*dSj  ^ 

Combining  (11-6)  and  (II-/)  yields 

if  - -t^F  /at) 
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Since  OCTFkWW)  > 0(rPRop^O(j^)~loOCr^.  Therefore,  ~ 


T.  - 3'. 

r V 


For  opposed,  parallel  disks  like  S(ol  and  S(x)  Sparrow  and  Cess ^ g|ve 


S(©}-*  S(X} 


llC  r„%  >f ) - f(Q%r,Vx‘)  - 4r;r„l  j /C2rf) 


ai-9) 


If  the  terms  in  (II-9)  are  non-dimensionalized  with  the  particle  diameter  D, 
the  form  of  the  equation  does  not  change.  Therefore,  (11-9)  can  be  employed 
in  nondimensional  form  (denoted  by  an  overbar).  For  a circle  with  center 
at  (r,xM  = (O,  R ) 


r*t  =1^ 


(IT-101 


Thus, 

r - x'(i- *<) 


so  that 

t C - t i\-L) 


(11-11) 


(II-  1 Z 1 


and 

rxl  = + 

Substitution  of  (11-12)  and  (11-13)  into  (II - 9)  yields  (with  difficulty) 


Fsbwsw  ' 1 . osK<  i-L) 

Therefore,  with  ( II  — 8) 

The  number  of  oxidizer  particles  on  the  burning  surface  vSp 
with  h and  D is^^) 

\c,/CTTDmV]  S^Ai 

Since  the  area  of  the  intersection  of  one  of  these  particles  with  Sp  is 
= TTr1  * -rr^ 


the  planar  surface  area  of  these  particles  at  the  burning  surface 
(x-o)  is 


d\  --  c,  l dl 


(11-14) 


(11-15) 


(IT-16) 


(II-  17) 


(11-18) 
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Therefore,  the  radiant  energy  entering  these  particles  is 


a*  9.  - a! 

(11-19) 

The  fraction  of  this  energy  deposited  at(&-D)£X$0  is 

d5a  - dF  /a 

0 S(p\-*dS3  0 

(11-20) 

With  1-15,  1-19,  and  1-20 

-3^  ■£-»  4^  dL  4x/d 

(11-21  ) 

Since  T = ^ y Sp  and  ^ - ciX/clx. 

<l*$  - 3D  dl* 

(11-22) 

To  deposit  energy  at  depth  -uX^O  # Therefore,  the 

energy  deposited  at  x by  all  oxidizer  particles  with  diameter  D is 
t CkkX* 

d§  = 3C>"  [ di1 

o 

(11-23) 

Consequently, 

d<$  - -U^s<o 

= 0 x 

(11-24) 

It  has  been  demonstrated  elsewhere^"^  that  the  particle  size 
distribution  in  an  oxidizer  mode  is  nearly  log  normal.  Thereto re^^, 

V pk2w 

|l  = i 

(11-25) 

Since  only  those  particles  with  IM  can  deposite  energy  at  x 

Ik*  1*1 

Note  that  unsteadiness  in  is  associated  with  that  of  J^, 

(11-26) 

■ 
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Figure  2 sketches  the  general  phenomena  involved  in  the  combustion 
process.  The  process  is  divided  into  two  major  parts:  the  nonreactive 
condensed  phase  (x  < o ) and  the  reactive  region  (x  > o).  The  latter  is 
depicted  as  containing  two  subregions- -a  gas  phase  reaction  zone  and  a 
condensed  phase  reaction  zone,  As  demonstrated  by  Novozhilovf 1 3 ' and 
Summerfield,  et.  al.  , the  characteristic  time  of  the  reactive  zone  has 
lesser  order  of  magnitude  than  that  of  the  nonreactive  condensed  phase. 
Accordingly,  when  the  order  of  magnitude  of  the  characteristic  time  of  any 
transient  process  is  greater  than  that  of  the  reactive  zone,  the  reactive 
z°ne  will  respond  in  quasi-steady  fashion.  That  is,  "inertia"  will  be  con- 
fined to  the  nonreactive  condensed  phase.  The  upshot  of  this  is  that  the 
reactive  region  may  be  described  by  appropriate  steady-state  relations. 

Moreover,  if  the  radiative  extinction  length  of  the  gaseous  products 
and  the  condensed  phase  is  large  compared  with  the  gas  phase  and  condensed 
phase  reactive  zone  thicknesses  respectively  , radiation  will  not  interact 
significantly  with  the  reactive  zone.  This  seems  to  be  a very  plausible 
assumption  for  reduced  and  minimum  smoke  propellants.  This  means  that 
for  a specified  formulation  the  reactive  zones  mass  flux  is  defined  solely 
by  its  boundary  conditions  (Ts  f)  and  environment  (p).  In  short, 

= rft\*r  ( formulation  ,1^4,  f)  (11-27) 

where  ( ) denotes  steady-state  conditions.  If  the  reactive  zone  is  quasi- 
steady 

'■  m-281 


and 


pu.  - 


-pu. 


(11-291 


Thus,  for  a quasi-steady  reactive  zone  the  instantaneous  burning  rate  is, 
for  fixed  formulation,  functionally  related  as 

^ = t)/o 


(11-301 


Note  that  if  radiation  does  not  interact  with  the  reactive  zone  that  this 
relation  holds  irregardless  of  the  radiant  flux  level. 


It  the  extinction  length  is  small  compared  to  these  dimensions,  the  geometrii 
scale  of  the  motor's  cavity  will  not  significantly  influence  the  radiant  flux 
incident  on  the  reactive  zone,  and  the  incident  radiative  flux  becomes  another 
independent  variable  for  mg.  Reference  20  has  attacked  this  problem. 
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Figure  IT-2.  Sketch  of  Regions  in  Combustion  Process 
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In  principle,  one  may  obtain  both  U-°(f;Te')  andT^°0(^T0)  from 
steady-state  experiments.  In  the  absence  of  radiation  an  energy  balance 
(steady- state)  for  the  x = 0 to  x = oo  region  yields 


* C - <C  (C -t.) 


(11-31) 


With  T^CfcT^known,  Eq.  (II- 31 ) givesT,  ■£).  Therefore,  u.“(y  T6\ and  T°  ( b T \ 
can  be  rewritten  as  J ° 4,0 


Clearly,  if^#{*  are  specified, is  also  defined  and  u!0  is  determinable 
from  Eq.  (11-30);  that  is,  from  the  reactive  zone  solution.  Therefore 
Eq.  (11-32)  is  quasi- steady. 

Cnder  the  optical  thin  assumption1,  the  reactive  zone  is 
independent  of  radiation.  Therefore,  the  connections  must  become 

“U*aV  u?u,n  , 

Finally  since  Eq.  (11-30)  holds 

■ u-u.o,  T^,(v,o  = \(^ 


becomes 


If  there  is  radiation,  the  energy  balance^steady- state) 


#*"-  w-T'C  - u + i/pcu:)] 


Therefore,  since  the  connections  have  been  demonstrated  to  be 
associated  with  the  reactive  zone  and  that  zone  is  not  effected  by 
radiation,  it  makes  no  difference  whether  f®  or  {®  is  employed.  With 
equivalence  of  C and  P Eqs.  (11-31)  and  (11-35)  show  that  radiation 
may  be  accounted  for  in  steady- state  by  employing  the  equivalent 
initial  temperature 


X*  ■ \ + 1/ 


In  other  words,  in  the  steady-  state  radiant  energy  deposition  and 
initial  temperature  change  produce  equivalent  effects  and  Eq.  (11-36) 
demonstrates  the  equivalent  temperature  for  specified  T0  and  J. 


(11-32) 


(11-33) 


(11-34) 


(11-35) 


(11-36) 
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The  equation  governing  the  non  steady  thermal  tield  in  the  non- 


reactive condensed  phase  (x  < o)  is 


“AT  o \ \ » » . x 

p*-s»  +f 


v It 


(11-37) 


The  boundary  conditions  are 


T(-ooyt^  * T0  5 I (<>,*}  - X, 


(11-38) 


and  the  connections  are 


'a,  ' t) 

where  ■f  1 C^T  fa*)  Introducing  dimensionless  variables 


(11-39) 


*»o 


~0--\Llo!t  £=U Jr./*,  £*(U.°)  t/x,  7 = t/V*  s ^ a f/f®,  f s f/-f’  (11-40) 


yields 

^B  *b©  ^ ^ , yc 

it  3 * +1 

with  boundary  conditions 


HI-41 ) 


* (T6-T*y(VT*)=  -OL  , --  15- 


(11-42) 


and  connections 


(11-43) 


Assume  that  small,  periodic  pressure  and  radiation  changes  are 
imposed.  That  is  assume 


v.'fl 


a tv  c* 


^ X e , ~4>  --  l + e. 


(11-44) 
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(11-45) 


where  0(1^*^  < 1.  Then,  in  response  to  these  tluctuations 

©Cs/t^  = ©°(s)  v ^ 1- 

'Dcc'i  ® i * <- 


Substitution  ot  these  iunctions  into  Eq.  11-41  and  neglect  of  higher  order  terms 
yields 


^ 1>V  "dB  ^ a /^e  is  oG  ^ „<■«*. \.**  d«*\ 

>Ut  " rsl  + M 4^n"e  fhdVdV 


J (11-46) 


Since  this  equation  must  hold  lor  >j,  - r 03  vji  =c  and 
dV  de°  vo  _ 

with  boundary  condition 

GC (-oo)  * -0L  , (d° i'-)  - b~° 

and 

^ „ is* 

with  boundary  conditions 


=o 


and 


*®.*  ^ ^ ^ VPl/t 
to  r ’ % *% 


with  boundary  conditions 


#• 

i4> 


(11-47) 


(11-48) 


(11-49) 


(11-50) 


(11-51  ) 


(11-52) 


Thus,  in  the  linear  context  the  problem  has  been  decomposed  into 
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Three  separate  problems:  (a)  steady-state  burning  with  radiant  deposition, 

(b)  burning  with  constant  radiant  deposition  and  nonsteady  pressure,  and 

(c)  burning  with  constant  pressure  and  nonsteady  radiant  deposition.  These 
will  be  considered  in  turn. 

Equation  11-47  can  be  rewritten  as  the  first  order  differential 

equation 


( ^°\  _ / 4b°\  _ y.* 

The  solution  of  this  equation,  noting  that  d6°/cJ§  - I at  g -o 

di  = e ~ e 

where 


> is 


(11-53) 


(11-54) 


= ^ r« e 


(11-55) 


When  there  is  no  radiation  ( ^ :0  ),  and  the  Michelson  distribution 

is  recovered. 


Substitution  of  (11-54)  into  (11-49)  yields 

* le.. 


at  " 'Sg'* 

a i/'t 


* n i 

TV  - ft  l)  ♦ ft  U Q, 

O % ,,yj  tV| 


Let  ( -(  ) ft  , then  (11-56)  becomes 

a p4c\  . a,  Sa  ( ^ 

-—"I  _ _U^Q  - e.  u -eu 

d%  "l  ,v>  "I 


d^ 

Pi  + 0P2  where  , ©p(  > and  6 are  solutions  of 


Take  © « ©H*  0_  + ©, 


d©,  d©, 


•:  - ^ -uSl.  =o 


,IA 


& ' If  * - - «' 


'"l 


- "i".  - i.V'B..  » £ eS  G 


This  form  puts  radiant  effects  into  ©f 


(11-56) 


(11-57) 


(I  I - 5 8 l 


(11-59) 


(11-60) 


Pi  • 
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respectively.  The  general  solution  of  11-58  is(21) 

©u  - A e tBe 


(Il-bl 


where  are  r°ots  of  the  characteristic  equation 


5, 


(ri-62) 


Since  must  be  bounded  at  § * -co  , B-O  and 

u 


9 * Ae 


(II- 63) 


where  ^ is  complex.  Algebraic  manipulation  and  use  of  the 

quadratic  formula  yields 


& - V -\Vs 


(11-64) 


(11-65) 


Th®2i?lution  (H-59)  can  be  found  with  the  method  of  undetermined 
coefficients  . Let  ©pi=C«  where  C is  the  undetermined  coefficient. 
Substitution  into  (11-59)  yields 


C = L 11  /t 


(11-66) 


Whence 


0p  * iVe  /t 


(11-67) 


Sin 


i e 


^ (note  that  H(o)>  = Q to)  =o) 


* 


By  definition  f /-T  = [(be/aSVC^&Va  6)1  . From  11-54  and  -55 

t - I . Thus, 


(11-75) 


$ 1+<^ 

Differentiating  11-74  (employ  11-73  and  Liebnitz's  rule  to  find 
(dH/d  S)^ _o  ~ O and  substituting  into  11-76  yields 


(11-76) 


A3.  +i*^^ 


(11-77) 


0_  »,F°».Ste!d-»'St«te  bu/ning  without  crossflow  the  functional  dependencies 

^ e ) an<3  ^ ^ can  ke  established  experimentally.  Therefore 

the  functions  9 


r^C'^/vr.*) 


IL 


' C Jo  ~ '*  -)  . £ (u  / iil/  p ) 


(11-78) 


can  also  be  established.  To  utilize  the  Z-N  method  one  needs  derivatives 
with  respect  to  and  because  of  the  universality  of  the  functions 

t)  and  L'.  ^.  T ) for  both  steady  and  nonsteady  conditions.  Since 


f~fC>  n^x)  / ^.T~f 
r «u\  c-^  . , w:  t. 

conversion  into  the  "universal  form"  requires  the  partial  derivatives 
U *-•+ , ?T  and  U-U.*  If  “>►)  * . Equation  11-35  gives  {Vuf  T T*\ 

Differentiation  of  11-35  then  yields  (with  11-78) 

' X ) C'-  4°/  = fc.  + r*-  i 

r 


(11-79) 


(11-80) 


' .0  \ •£- 
^ <3  • T , ■ ^ J - ^ t 44 

T* 
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The  Jacobian 


X - Aj(x.°  ) / '-(T^  - v?  k. 


(11-82) 


Since  and  B-- ) [dimensionless  form  of  andT^i't*) 

t;  cr4^] 


13,  * (dD/o<^  $ -13,  'c^  ^ 


A -c'c^/^^4,  ^-C  ^1. 

With  the  definitions  for  ^ ^ 'v]1  ^ 
(dtt/uvj)^  * 0)/yj  ) '■ 

Cci9-/ >-Y[^  --  ^ 

c>  >/  0,  ‘ (r:  r; )] 1 Col; / b U\ 

t r 
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ai-83) 


(11-84) 


(II-8S) 


(11-86) 


Therefore,  neglecting  second  order  terms 

+•  fe, 

"&,*  = r*<f,  /(  f U’+yu.*)^,  /(  k\r*-i) 

With  11-75,  -77,  and  -86  there  are  four  equations  and  five  unknowns 
A ).  Consequently,  the  ratio  T\  / Tq , can  be  solved  for  in  terms 
of  quantities  that  can  be  computed  from  uf  (>,T0*)  and  7^  ( V.  C*)  data. 
Since  i),  /-q,  is  precisely  the  pressure  coupled  response  function,  this  is 
the  result  desired. 

Examination  of  11-75,  -77,  and  -85  show  that  none  of  these  equations 
depends  explicitly  upon  the  radiant  energy  deposition.  Consequently,  it  must 
be  concluded  that  the  effects  of  steady  radiant  energy  deposition  in  the  non- 
reactive solid  can  always  be  accounted  for  by  employing  the  radiation 
augmented  initial  temperature  T0  ‘ j That  is  the  pressure  coupled  response 
at  that  at  >,T0j^-o.  This  situation  was  deduced  on  physical 

grounds  for  the  special  case  £ This  work  shows  this  to  be  the 

case  in  general. 
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APPENDIX  A 


NOMENCLATURE 


Latin  Symbols 


constants 


e,  exp 


specific  heat  or  acoustic  speed 
constant  or  function 
diameter 

weight  mean  diameter 
denotes  the  error  function 
denotes  the  exponential 

temperature  gradient  at  x=o,  (dT/dX)x_Q 
direct  interchange  factor  for  radiant  transfer 

f e.  * d( 

O 

see  Figure  II-  1 

o o 

tf-l1  or  an  index 
denotes  an  integral  (see  1-13  to  -16) 
radiant  heat  flux 
radiant  heat  flow 
an  index, 

function  defined  by  11-78 

extinction  depth  for  thermal  radiation 

interaction  parameter  (see  1-3)  or  mass  flux 
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r 


M 

N 

n 

P 

r 

R 

R 

P 

Rv 

S 

t 

T 

t* 

T0 

ri"' 

W 

X 

X 


(Cont.  ) 

number  of  oxidizer  modes  in  propellant 
number  of  oxidizer  particles  on  surface  Sp 
pressure  exponent 
pressure 

radius,  burning  rate 
function  defined  by  11-78 

see  F igure  II-  1 

pressure  coupled  response  function 
velocity  coupled  response  function 
surface  area 
time 

temperature 

T+  J* 

initial  propellant  temperature 

burning  surface  temperature  or  burning  rate 

<*-*  )/(<r  fT  ) or  burning  rate 

mass  fraction 

InD  or  spatial  coordinate  normal  to  the  mean  burning  surface 
see  Figure  II-  1 
roots  of  11-62 


43 


APPENDIX  A 


(Cont. 


Special  Symbols 


(To  -T*  )/(Tj-i;*) 


mass  fraction  of  oxidizer 


K /(U.0)2 

defined  by  II-  82 

dummy  variable  of  integration 

dummy  variable  of  integration  or  temperature  gradient  ratio  f/f 
density 

function  defined  by  11-7  8 
radiant  energy  deposition  per  unit  depth 
nondimensional  temperature 
thermal  diffusivity 
burning  rate  ratio  A.° 
non-dimensional  distance  w / v. 
non-dimensional  time  (u?)vt  IK 
non-dimensional  pressure  ^/>° 

non-dimensional  radiant  energy  deposition  $ / $ 
volume  fraction  of  oxidizer 

non-dimensional  radiant  energy  deposition,  $*(K/0/[pcAt -O] 
function  defined  by  11-78 
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(Cont.  ) 


*• 


( ) 


o 

( ) 

\ 


non-dimensional  burning  surface  temperature, 

^ or  standard  deviation 

denotes  a spatially  mean  value  or  a value  non-dimensionalized 
by  the  oxidizer  particle  diameter 

denotes  a pseudo-propellant  property  or  conditions  with  J^=o. 
denotes  steady-state  condition 


Subscripts 

b denotes  binder 

d denotes  diameter  dependence 

i denotes  ith  value  or  imaginary  part  of  complex  number 

k denotes  kth  value 

ox  denotes  oxidizer 

I denotes  a perturbation  quantity 

denotes  the  portion  of  the  perturbation  quantity  belonging 
to  the  pressure  perturbation 

Ity  denotes  the  portion  of  the  perturbation  quantity  belonging 

to  the  radiant  flux  perturbation. 

t denotes  total 

p denotes  quantity  based  on  planar  area  Sp 
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APPENDIX  A 
(Cont.  ) 


o 

r 

s 

S 


denotes 

denotes 

denotes 

denotes 


Js=o  condition 
the  reactive  portion 
conditions  at  burning  surface 
the  area  the  function  pertains  to 
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APPENDIX  B 


FORTRAN  IV  CODE  FOR 
EXTRACTING  PSEUDO-PROPELLANT 
RATES  FOR  SPECIFIED  M 


□ I *E  NS  ION  M(8  ) » HE  AO  (201  tPf  5)  ,RMA<8  ),  I PC  INK  16) 
CQ^MGN/MADER/lI.JJ.At  F A ( 8 .SQ)  iR  (50  *_5  ) . Ml  50  ) , RC  t 5 Q ) , AL  FA  113  3 )_»_ 


COMMON  7PLK/dM0t 6) .SDltil 


K.i(8f5)fNC(50)tKfEpp.  (50)  , ERNC50) 

A1  (8)«-AJlZ51»AK(8J  tUllfllt. 


00000020 

00000010 

00000043 

QGQQ0D5-Q 


XMAX( 50) ,XMIN(50),RH0T(  50) , RHOX , RHOB , AP I , SEEJ,  00000060 

C ( 5 J ) , r<  A ( 1 7 ) « Nl  a AN  ( 6 ) , DI_A  t S IG  ^SPLN  ( 8 ) . h'MDLN  I Q.L , IQ. QCQOOQLQ 


COMMON  /MAN/  E 1(8, 15)  , F2 ( 8 , 15 ) , F 3 ( 8 , 1 5 ) , E4 ( 8 , 1 5) , EL EMI ( ? 1 , 8 , 1 5 > , 000000  8 0 

1 ELC"2(21>8.  15).  EL  EM  3(2  1,8,  15)  t EL  EH4(  2 1,  8 1 15)  .XC16) ,12  000Q0Q2D 

REAL  N,NM,NC  00000100 

C R ( J «K  1 EXPERIMENTAL  BURN  RATES  FOR  THE  JTH  FORMULATION  ANC._Pm QQ  QJQ113 

C N ( J ) EXPERIMENTAL  PRESSURE  EXPONENT  FOR  THE  JTH  FORMULATION  AND  ^0000120 

_£ PJJJ 00000130 

00000140 
QC0QQ15D 


C PC(J)  CALCULATE 0 BURN  RATE  FOR  THE  JTH  FORMULATION  ANO  P(K) 
C NC(J)  CALCULATED  EXPONENT  FOR  THE  JTH  FPPMULAT1QN  AND  PULI 


CKRi  STANDARD  DEVIATION  OF  RC  FROM  k FOR  P(K) 

ERR2 STANDARD  0EV1ATIQN  OF  NC  FROM  N FCR__P(K_L 


00000160 

QCQQ0.1I3. 


R M ( I , K ) MODAL  BURNING  RATE  FOR  ITH  MODE  AND  P(K) 
NM(t)  PRESS UKh  EXPONENT  FOR  THE  ITH  MODE  ANQ  P(l) 


00000180 

00000109 


C ALF  A{ I , J ) MASS  FRACTICN  Of  ITH  MODE  tJXIDl ZER  IN  JTH  FORMULATION  00000200 

C ALFATCJ)  TOTAL  OXIDIZER  CONTENT  QF  JTH  FORMULATION QQQQQ2JJ 


C II  NUMBER  OF  OX  MOOES  < 8 

C JJ njmhfr  u f formulations  <,50 


0C0002  2 0 
300002.10 


C KK  NUMBER  PRESSURES  RATE  DATA  AVAILABLE  < 5 

C ERR ( J ) ERROR  BETWEEN  RIJtK)  AND  PC(J)  AT  P(K) 


00000240 

000QQ250 


C ERN  C J ) 


ERROR  BETWEEN  NIJ)  AND  NC(J)  FOR  P « 1 ) 


00000260 

OQOOQ21Q 


C THEORY  ASSUMES  THAT 

_£ K (J  ) = Su.M  f ALFA  ( I , J ) «RM(  I ) )/  SUM  t ALFA  ( I , J)  ) 


00000280 

00Q0Q23.Q. 


N ( J ) = SJM ( ALF A ( I , J)*RM< I ) *NM ( I ) )/ SUM ( ALFA ( I,  J)*RM.(  I ) ) 


00000300 

00000110- 


KEN  HERR EN  MARCH  27,  1978 


00000320 

030  QitflQ 


INPUT  J0D  HEADING  AND  'PRINT  OUT  THAT  HEADING 


00000340 

QQQQQ11Q 


.Ml  *6 

APIO. 1415V 


00000360 

QDQQ.QULQ 


RHOB=0. 91 
K HCX  = 1 » 9 5 


OOOOC3BO 

QQQQQ12Q 


1 READ(5,10, FND=210)  HEAD 
10  FORMAT ( 20A  4 ) 


000004D0 
OOP  0Q4.10 


WRITE(6tll)  HEAD  00000420 

11  FORMAT! lhl  ,2QVt,//2H  »? OX.«OX  MASS  FRACTION  OAT  A*  , 2 IX , 2 H*»  , 2 IX  0000Q41Q 

l, • EXPONENT, RATF  DAT  A • , 2 3X , 1H* , / / MX , • MODE  1,,2X,,M0DE  2 f , 2X , • MODE 3 • 000 00440 
2 , 2 X « * mtj  jE  4*  , 2 X « * MODE  5».2X.«MC0E  6«,2X,'M0PE  7«,2X,«MQDF  8«»3X,  00000450 


3 #N (P 1 ) #, 3X,»R( PI  ) f,  MX,»R(P2) * , 3X, »R(P 3) * ,3X, »R(P4)'  ,3X, «R(P5)f ,/) 


00000460 

QDQQ1VZQ. 


INPUT  '‘A  T A 


00000480 

000004.90 


RC AJ( 5,20)  I I , J J , KK , ( P(K) , K= 1 » KK ) 
20  F n R M AT ( 312  .SF5.0) 


00000500 
OOP 00510 


DO  60  J=l, JJ 

R E AD  ( 5 , 30  ) ( ALFA(  I , J)  , 1 = 1,8)  ,N  ( J)  t (R  ( J,K),K=UKK)_ 

3D  FORMAT! 14F  5.4) 


00000520 

Q3QQQL3Q. 


00000540 

QQQQQ5-5Q 


OUTPUT  INPUT  DATA 


0000056  1 

QPQ935J3 


WRITE(6, 40  ) IALFA! ), J), 1*1, II  ) 
40  FORMAT  C 1H  ,5F8.4) 


WRITE (6, 50  ) N(JI,  (k(J,K),K=l,KK) 


0O0005B0 

00000600 
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U o KJ 


25 

FORMAT ( RF7.3) 

00000640 

R E AD  ( 5 , 2 7 ) (SJ(I).l  = L.II) 

2 7 

FORM  AT ( 8F  5 • 3 ) 

00000660 

WRITE (6, 2 6 ) ( WMU(  I ) , SO  II) , I = 1 . I I) 

26 

FORMAT( lXf 2F10.3 ) 

00000680 

RE  AO ( 5 * 29 ) (RA( I ) , 1=  I,  13) 

00000690 

2^ 

FORMAT (3E 1 2 #4  ) 

00000700 

READ ( 5, 31 ) I I POINT  ( ‘IN  ) * NN-2  *11)  ,12  t I JKLM 

00000710 

31 

F OPM  AT ( 1 2 I 2 ) 

00000720 

WR  I T £ ( 6 f 32  ) R A ( J 2 ♦ 2 ) 

00000721 

32 

FORMAT ( 1 HO , * THE  CURRENT  N = *,F8.4) 

0000072 2 

C 

00000730 

C RESET  KPRINT  FCR  NMII)  DETERMINATION 

00000740 

C 

- 00000750 

KPRI NT= 1 

00000760 

C 

00000770 

C COMPUTE  ALFAT ( J) 

00000780 

C 

00000790 

00  6 7 J=1,JJ 

00000800 

ALFAT ( J ) =0 .0 

00000810 

DO  63  1=1,11 

00000820 

ALFAT < J )= ALFAT ( J ) +AL  F A ( I , J ) 

00000830 

R HOT ( J) sPHOT ( J ) ♦ ALF A ( I , J)/RHOX 

00000840 

63 

CONTINUE 

00000850 

KHOHJ  ) = 1 • / ( RHCT  ( J )M  1. -ALFAT  ( J ) ) /RHOB) 

00000860 

67 

CONTINUE 

00000870 

C 

00000880 

C FIND  RPUvKI  FCR  EACH  PIK) 

00000890 

C 

00000900 

00  230  K = 1 , K K 

00000910 

C 

00000920 

C SET 

INITIAL  R M ( I , K ) 

00000930 

C 

00000940 

WR  I T F ( 6 , 68  ) 

00000950 

68 

FORMAT  ( 1 HO,  • MODAL  INTEGRALS  AREf /I IX,  *XMAX*  ,11X,*XMIN*, 

00000960 

1 8X, 'DIAMETER* , 8X,* SIGMA* ,8X, 'INTEGRAL* ) 

00000970 

00  70  I * It  II 

00000930 

WMOLNl I )=ALOG(WMD( I ) ) 

00000990 

SULN ( I ) =ALCG( S0( I ) 1 

00001000 

X M A X ( I ) =WMQLN ( I )+3.*SDLN( I ) 

00001010 

X M I N ( I ) =WMOLN ( I )-3 • *SDLN( I ) 

00001020 

XY1=X*AXII>  _ 

00001030 

X Y 2- X Ml N III  p 

00001040 

3I*V»3LNm 

00001050 

s I c=  SDL  M ( n 

00001060 

AI  { I ) - S I M P ( X Y 1 ,XY2,A2,Nl) 

00001070 

00001080 

69 

00001090 

R MA ( I ) = P ( 1 , K ) 

oooonoo 

IKK  .EQ.  1)  AM  I )=N(  1) 

00001110 

70 

CONTINUE 

OCO  0 1120 

C 

00001130 

c 

COMPUTE  THE  FORMULATION  INTEGRALS 

00001140 

c 

00001150 

X S r A R T = a M D L N ( i ) +0.U  3. * SDLNl  l) 

00001160 

XLNJ  = W MO  L N ( 8 ) -0 • 1 - 3 • * S DL  N ( 3 ) 

00001 i 70 

XDEL  = ( X E NO—  X S T A P T ) / 3 1 • 

00001)30 
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X(IL)=X$TAkT  + FLQAT(  IPQ1NTC  IL  ) )*XDFL 
vt  f mi  it- 


IL T= I 2+  1 

1 Q-=  1 LT.tJ 

X C I L T > = Xf  M) 


838  FORMAT  ( 1M0  ,P(  ?X,F  1?.  5M 

L*  JJ 

SUMaO. 


iCj  0 I = L t 12 

i = y ( T ♦ 1 ) 


IF(ALFA<K*=\*.J).EO.O.  ) GO  TO  2100 


[)  I A = w ML)L  N ( < E N ) 


niKf-Nt  I )=—  SIMfMYlfY2tAl»M) 
F?(KL'l\.  I )=-SI  4P(Y1, 


tw 


00001220 
0 *)  A 1 


OC001240 

OGOQ125J 

00001260 


00001230 
JC  0012210 
00001300 


000  01 3 ? 0 

' ITT 


"00013^0 


00001360 


00)01330 
00 001300 _ 


IMA71.1F.50.0)  AZ  = EX'P  (-0. 5*  A 7 1 ) 
I F ( 3 Z 1 • L L « 5 3 « 0 ) B Z = tL 


f 3 ( KENt I > = BZ-AZ 

4 ( K F M • 1 ) = *Z»FXP<  (RAC  1G1-3. )»X 


A7*EXPURA(IG)-3#)*XtI*l)) 

_Z= X ( I »1 )-X ( I ) 

T = ( X ( I )-V  v :)  L N ( KEN)  )/Z 
FLE  MIC  J tKRJ 


*ALFA(KENt J ) / 2 . 506627 
f if  ' V ( J fKrM»  1 ) = ( — 1 • # T # E 1 ( KENt  H/SDlNf  KEN 


UFACKCNt  J )/2. 506627 

Li-  ^7..(.J.tii£N-».U-=  ( ( 1 .+T-(  (i~  .MIGl-3._i»SDLN(KEfl>  )/Z)*E2HlEI 
/$  >L:\'(  KKN)  4-  F4CKEMtI  ) /Z ) *ALF  A ( KEN  ,J  ) /2 . 506627 
F L 1 m 4 ( J . K F NJ , I ) = ( - 1 • * ( T - ( (RACIG)  — 3#l  *$DL  K!lK£N 


E2  ( <V*  # 1 ) / V)L  N(  XF  N)  — E '♦  ( KFN * I ) / Z ) *ALFA  ( KFN * J )/ 2. 50662  7 

or  TO  3 50  0 I 


iMagum 


L L T I \nr 
r L!  Ml(  J,KL-Nf  1 ) 
LLV--.2CJfK.tNf  l ) 
L Lf  MU  J • K ‘ I , I ) 


±_LLzJ.__ 

.1)0. 

. I ) =0. 


r S'  « % 


~SirAi'X2 


(,  OMTINUE 

AJ(.i  i = SUvV2.^0662  f 


00001440 
O 


00001460 


00001430 
03^  1A9:Q_ 
00001500 


00001520 
Z)»000 


00001540 

I)  QQQ-Q  15.53 

00001560 


0 3 0 J I 5 3 J 
00001590 
00001600 
00  OJ  liU 
00001630 

oooo  mo 


0C0O  16  30 


16  49 


1000  CONTINUE 


WRITE  (6.  12!45)  ( I El  (KCNTt  ICNT  ) , ICNT=i  t 1?)  , KC  NT=  1 , I I ) 
RpiT£!6,1235)  ((E2(KCNTt  ICNT) , ICNT=1 , 12) t KCNT= 1 » I I ) 


IwRI  T C ( o , 1235)  ( ( E3 ( KCNT  , ICNT)  • ICNT  = 1 ,12)  , KCNT=1 , I l ) 
WRITER,  1235)  ( (EX  KCNT»  ICNT)  , ICNT=1 .12)  ♦ KCNT= 1 ♦ I I ) 


1235  Ff,P M AT  ( 1 H ,IIF8.4) 

C 


1234  CONTINUE 
C PINO  UPT  I M A L R;M  I t K) 


00001780 

0000179 


00001810 

00001820 


00001830 

00001840 


00001850 

00001860 


00001870 

00001880 


00001890 


CALL  PAT  SHI  AN, ERR 2, I I f5.0f 0.001 


DC.  95  I = l,  II 
95  NM  ( I ) = AN  ( I ) 


C 

C OUTPUT  COMPUTED  RESULTS  THIS  p{K) 


100  WRITE(6,UC)  HE A'O 


110  F 0 P * 1 A T ( 1HL  , 2 0 A 4 ) 

WR  I T E ( 6 , 1 1 5 ) PIK) , FRR1, EPP2 


115  FORMAT ( 1H0 , 3HP  =,P6.1,4H  PSI,5X 
l NT  STD  FRP =• , IPC1  1.4///2H  *,16X 


22H**  tl3X,  ’ CCMP -JTCD  ^00 A L EX^ONE 
3,  • RM  3*  , 5 X , * R M 4 • ,5X,  16*  ,5Xf 


4'  Ut*2  9 ,5X  , • \.\3  • ,5X  , •NM4f  , 5X  # »<;M5 

WR1TF  < 6, 12C)  (RKI  I ,K)  , I =1  , I I ) 


120  FORMAT ( in  , RFR.4) 

I F ( KPR I NT  • F 9 • 1)  WRITE (6, 130) 


130  FORMAT! 1H>,64X,  HP  8. 
WRITE (6, 140) 


Oi*  ^AT  ( 2H  ')  v t I ,)x  , * C OMP  AP  I SO  J THE 
1 • R C * » 5X , • ERR • ,6X, •N* , 7X, •foC* , 5X 


DO  1 80  J=1 , JJ 

wRITE(6,15  0)  P(  JtO  yRC(  J ) , E H 

150  FCPMA1SIM  j -i  F f,  • 4 ) 

I PI KPR I NT  .FO.  L)  w P I T P ( 6 t 16 


FORMAT  ( 1H4-  , 24X.3F8.4) 
180  CONTINUE 


KPPINT=1 
200  CONTINUE 


CiU  TO  1 
210  STOP 


00002010 

» 5 0,-1 , f‘CUFf  2 ) 0 0002020 


00002030 

00002040 


00002050 

00002060 


00002070 

00002090 


00032090 

00002100 


, * R A T r STD  EPP  = » f 1PP1  1.4, 5X, 1 EXPONE00002 1 1 0 
, • COMPUTED  NODAL  P ATE S, I N/SEC • 9 18X , 00002 l 20 


NTS'  »2lX,lH*,//4X, •RMl,,5X, 1 K K 2 1 , 5X00002130 
6 1 , 5 X , • R N7 • 9 5 X , • P M 8 * , 5X f ■ NM1 1 , 5X , 00  3 02140 


• , 5X,  • N*6»  , 6X9  'M7*  , 5X9  • NM8f  , / /)  00002150 

00002160 


00002170 

( ( n , i = i,  m 00002180 


00002190 
00°02?  00 


rR  Y/ EX  PER  I VENT • , 9 X , 1 H*/ / 5 X , • 2 • , 7X,  00002213 
, • ER  N • ,/)  00002223 


00002230 

R!  J) 03302243 

00002250 

0)  N(J),KCIJ),ERMJ)  0CO02  260 


00032279 

00002280 


003 02293 
00002309 


00002319 

00002323 


o olo  nln 


SUBROUTINE  MQl)EN?(  AN,  ERR2) 


EXTtPNAL  A1,A2,A3 
OIMENSION  WUMj) 


COMMON/ BLK/fc’MDC  8),  SOI  8)  , AI(R),AJ(25),AK(8),U(50), 

J XM  A X ( 5 C ) tXM!N(50>  t«HGT  (50)  , P HOX , R HOB , AP I , SEEJi 


REAL  N»NMyNC 


C( 5 3) , S A( 17) ,N1 ,BN( 8) , 0 1 A , S I G , SOLN ( 8 ) ,WMDLN(8 ) ,TQ 


QOOQQQ1 Q 


00000020 

00300033 


COMMON/ MAD  EK/  1 1 f J J , ALFA! 8 v 50 ) « R( 50«  5) «N ( 50 ) »RC ( 50) • ALFAT < 50)  00000040 

l ,RM(6,5I ,NC(50)  # K < ERR  ( 50  ) t ERA] (50)  00000050 


00000060 
■ QQQ0007Q- 


00000080 

00000090 


C 

c 


THIS  SUBROUTINE  FINOS  THE  STANDARD  DEVIATION  BETWEEN  NtJ)  AND  NC  I J ) 


N ( J ) 


FXPER  IMFN'TAL  EXPQNE  MT  FOR  JTH  FORMULA  T I ON  AT  P(i) 


C NC(J)  THEORETICAL  EXPONENT  FOR  JTH  FORMULATION 

C KM(I,K) BURNING  RATE  FOR  THE  ITH  MODE  AT  P(K) 


C 

C 

C 

C 


NMU)  EXPONENT  FOR  THE  ITH  MODE  AT  P(1J 

_A L FAU»J)  MASS  FRACTION  OF  ITH  OX  MODE  IN  JTH  FORMULATION 


ALFAT ( J ) 
I I 


TOTAL  mass  fraction 

NUMBER  CE  MODES  <8 


jF  QX  IN  JTH  FORMULATION 


c 

c 


JJ 

ERNU) 


NUMBER  OF  FORMULATIONS  < 50 
ERROR  BF.TaEEN  NC  J ) AND  NC(J) 


ERR? 


STANOAKD’EKROR  OF  ESTIMATE  FOR  DATA  SET  N AND  NC 


KEN.HEKREN  MARCH  2«,1978 


DO  15  I=l,° 


% 


00000100 

00000110 


00000120 

00000133 


00000140 

00000150 


00000160 

00000170 


00000180 

00000190 


00000200 

00000210 


00000220 
0 3 J 0 0 2 3 0 


00000240 

00000250 


00000260 

00000270 


15  BN l I ) = AN ( I ) 

5 ERR2O.0 


LO  30  TO  J=1 , JJ 
NC ( J ) = ).0 


V 


00  25  1=1,11 

ir  ( w m o ( i ) •r.r.o.)  go  to  23 


% 


00300280 

00000290 


00000300 

00000310 


00000320 

00000330 


7NCALC=  3. 
GO  in  ?<+ 


2 5 


Y X = X ' ' A X ( I ) 

YY  = X MINI  I ) 

!.)  I A = W M 3 L N ( J ) 
SIG=S  )LM(  I ) 


% 


/ K { I ) = S IMP  (VX  , YY  , A3,  91  ) 

TNC  A l C = A(  F A(  I f J)*AK  (I  ) /SDL  HI) 


24  CONTINUE 

25  I f ( JJ  =J  1 r ( JI+TNC/ 

NC ( J ) = NC ( J )*0.3  ;«94 
CALCULATE  tvpnrf.NT  USIN' 


*LC 


v 


ME  A SURE  3 PATE 


•!C(  J) 
CALC  l ) L M e 


: NC  I J l/'ALF  AT  ( J)  /R<  J,l) 

EXPONENT  USING  CALCULATED  RATE 


NC  I J)=\iC  ( J ) /ALFAT  ( J)  /KC I J ) 

CAJLCUL  AT  E _AHSNLUTE  ERF  UR 

lR: (J ) =N( J )-NC ( J I 
CALCUL  AT!'  ^FL  ATI  Vf  l RK  IK 


\ W ( J )= (NI JI-NCI J ) )/N( J ) 

3 0_J*_R  ' .’  = F ’ p 2 ♦ FPNI  J)*ERM  J) 
\ RR  ? = S J R T ( ER  R 2 / J J ) 


00000340 

00000350 


30000360 

00000370 


00000360 

00000390 


00000400 

00030410 


00000420 
000  004  3 0 


00000440 

00000450 


00000460 

00000470 


00000480 

00300490 


0000050.0 

00000513 


00300523 

00000530 


33000540 


SUr.PCUT  INl  MQJEP?  (RA,  £ p R 1 ) 

EXTFRNAL  A 1 , A2 , A3 

C G 'V  M t : N!  / .v.  A ;~j  [ R / I I , J J , ALfM  8. 50  ) , R ( 50,  5)  , M 50)  , RC  (5  0)  , ALFAT  ( 5Q| 


,&!*(8.5)  fNC<50)f  K,ERK!50>  f EPNI50) 

C 0 M.M  U N / h L K / a M tD  ( ?.  ) , S)(  8)  . AI  ( 3)  y AJ(  25)  ,AK(  8)  ,U(  50  ) , 

XMAXI 50) , X M I N ( 5 0 ) , RHOT (50) , RhOX , RHOB , AP I , SEE J t 

f.(  50)  , SA(  17)  ,N1  ,BN(  8 ) , 0 1 A , 5 I G , $ DLN  ( 8 ) ,WMDLNC8)  tTU 

c n iM  ( JM  / M A N/  Ei(8»l5)fE2(8fi5)fE3(8fl5)tE4l8fi5)fELEMl(21f8fl5)t 
FIE  M2 (21,3,15) t ELEM3 ( 2 1 , fl , 1 5 ) , EL EM4 ( 2 1 y 8 , 15  I , X ( 1 6 ) , I 2 


DIMENSION  KM  A ( 3)  ,RA(  17) 

THIS  SuPPCJTINE  COMPUTES  STANDARD  OEVIATtCiN  BETWEEN  Ri  Jt  K)  AND  RC  ( JJ 
FUR  P ( K ) 


WHERE 

R ( J , K ) 
k C ( J ) 


RMAC  l ) 
ALFA! I , J) 
ALFAT  l J ) 

M 

JJ 

E»P< J) 


EXPERIMENTAL  BURNING  RATE  FCR  JTH  FORMULATION  AT  P(K) 

T ^EfiRET  I CAL  BURNING  RATE  FCR  JTH  FORMULATION  AT  P(K) 

B JD  N I NG  "~R  AT  E FOR  ITH  MODE  OF  FORMUL  AT  I CN  AT  P(K) 

*ASS  FRACTION  OF  ITH  MODE  QX  IN  JTH  FORMULATION 

TOTAL  OX  MASS  FRACTION  IN  JTH  F CK  MULA  T I ON 

NUMBER  OF  OX  MOOES  < 0 

NUMBER  OF  FORMULA T I QMS  < 50 
ERROR  BETWEEN  R(J,K)  AND  DC(J)  AT  P (K ) 


TAN DAP D ERROR  OF  ESTIMATE  FOR  DATA  SET  R AND  RC 


KEN  HEKRF.N 


MARCH  27, 


ILT*I2»1 

DO  ? I * I , ILT 
5 S A ( I )=RA( I > 

10  ERR  1 = 0.00 
SUM1 =0 • 

DC  300  J= 1 , J J 
RC  ( J ) = 0 • 

S U M E = 0 • 

SUMP  = 0 • 

DO  330  1=1, 
SJMl =0. 
SUM 4=0 . 


SU  

SUNG  = 0* 

SUMO  e1» 

00  350  KEN  =1,11 


SLMA  * SUM A ♦ E LE  Ml ( J , KFA , 1 ) 
SL’B  = SUM  3 ♦ E L E M 2 ( J , K E N , I ) 
S LMC  - 5UVC  ♦ E LE  M3 ( J , KFN , I ) 


SUMD  = SUMO  ♦ ELEMM J,KEN, I) 

360 CUNT  I.NUE 

SUVE  = KA(  I )*SUMA  4-  R A ( 1*1  F* ' S U M B ♦ SUN 

SJMF  = RA(I)*SUMC  » R A ( I ♦ 1 ) ♦SUMO  ♦ SUN 

330  CONTINUE 

K C ( J I = SUM  c + < l.-ALFAT( J) )»SUMF/AJ(  J) 
1)3))  SUM  A , SUM'5 1 SUMCf  SU^i) , SUM  E t SUMF 

C 1000  FORMAT!  1 M >nx»M?X,F12.5)) 

C COMPUTE  ABSOLUTE  ERROR 

E K 3 ( J)=P(J,K)-PC(J) 

C COMPUTE  P EL  A T ! Vr  EPF  OF 

C FKRU)>(R(J,K)-MC(J))/R(J,K) 

EK-f  1*1  R«”l*FPR  ( J)*FRKj  jS 

C WRITE! 6 ,35)  CRR ( J ) , R ( J , K ) , PC ( J ) , ERR  1 

C 35  FL'LMAff  lh  ;4E12.5) 


RAC !♦!) vSUMB 
R A ( I 4-  X ) ♦SUMO 


00000020 

00000030 


00000040 

00000050 

00000060 

00000070 

00000080 

00000090 


00000100 

00000110 

00000120 

00000130 

00000140 

00000150 


00000160 

00000170 

00000180 

00000190 

00000200 

00000210 


00000220 
00000230 
00000240 
00000250 
00000260 
OOP  002  7 0 
00000280 
00000290 
00000300 
00000310 
00000320 
00000330 
00000340 
00000350 
00000360 
00000370 
00000330 
OOC  003  90 
00000400 
00000410 
00000420 
00 300430 
00000440 
00000450 


00000460 
00000470 
00000 430 
Q000Q490 
00000500 
30300510 
00000520 


00000530 
00000540 
00000550 
00000560 
OCO  00  5 70 
00000580 
0Q0Q059Q 
00000600 


53 


00000030 


r 


SUBROUTINE  PflTSHIPSI.SSI  .N.  DELS.  OLMIN.IT  II  M.  I PT.MEP  IT) HQQQOfli.Q  - 

c 

UI  ME  NS  ION  PS  I (.75)  ,PHI  (25)  * THT  ( 25  ) , XFI.G  ( 2 5 ) 

00000030 

00000040 

c 

PRINTFR  AN  O H I r»H  SPEED  CONSOLE  DFVICE  NUMBERS 

00000050 

c 

-CLQQD-QI16Q 

DATA  ALFA/1.02/ 

00000070 

c 

QQQQQQfiQ 

c 

FUNCTION  F IS  MINIMUM  IMPROVEMENT  REQUIRED  OVER  LAST  BASE-POINT 

00000090 

c 

!>!«■! 

c 

F ( SS  S) =SSS  - A3S( SSSJ*0.0001*CUT 

00000110 
-33300123 

c 

PSI  IS  THE  CURRENT  BASE  POINT 

00000130 

c 

THT  IS  THE  PREVIOUS  BASE  POINT 

— ilihlilill  liM 

c 

PHI  IS  THE  TRIAL  POINT 

00000150 

c 

Milil  III  IMiBi 

c 

IPT  = 1 PCK  DIAGNOSTIC  PRINTOUT 

00000170 

c 

D FTR  MINIMAL  PRINTOUT 

■ 1 II  Mil  ■ 

c 

-1  F CP  NO  PRINTOUT 

00000190 

c 

• 

c 

00000210 

c 

INITIALIZATION 

00000220 

c 

IT'aICE=1 

00000230 
QQQQQ23  5 

l)EL  = DELS 

00000240 

00  705  I®1 

00000260 

7)5 

X FLG ( I ) = 1 • 3 

,3,0000270 

ITEP=C 

000  002  BO 

C J T=  1 • 0 

c 

00000300 

c 

EVALUATE  AT  INITIAL  BASE  POINT 

c 

1 0 

CALL  MEM  IT (PSI •SSI) 

00000320 

00000330 

c 

00000340 

c 

'EXPLORE  AROUND  CURRENT  BA3FP0INT 

.20000230 

c 

9 0 

SSITST=F(SSI) 

00000360 

QQQ0Q370 

1)0 

S=  SS  I 

000003B0 

?i  P A T M = D 

Q0QQQ.29Q,  .. 

DO  l 01  1=1,25 

ODDOD4DO 

101 

PH  I ( I ) = PS  I ( I ) 

QQOQQ^.IQ 

I C AL  L= 1 

00000420 

IFCIPT.LT. 3>GQ  TO  15D 

00000430 

WRiTEIfrt  *5991  I TER 

00000440 

WRITE  (6, 60  0 ( PSI  ( J ) , J=1.N) 

m'lXN'iwm 

WR IT  E ( 6 , 60  l ) S , DLL 

00000460 

c 

c 

MAKE  EXPLORATCRV  MOVE 

00)00480 

c 

GO  TO  150 

00000500 

c 

c 

IS  THE  PRESENT  VALJE  LESS  THAN  THE  BASE  PCINT  VALUE 

00000520 

c 

. ■l.'l-i-rrY.  P-WW* 

160 

IF ( S.L  T. SS ITST)  GC  TO  200 

00000540 

c 

0Q0QQ55  0 

c 

CUT  STErj  SIZE 

00000560 

c 

— 1 II  II 1 1 IB 

DE  L = ).5*DEL 

00000580 

IF  IDE  L -01.01  MIN  ) GO  TC  100 

-QQQ-Q05-9Q 

IF  C 1 PT.Gf  . C)  a«I  TM  6 , 704) 

00000600 

55 


IHCUT.LT. 0.5)  GC  TO  702 

0000061 0 

c 

c 

ST  AWT  OVER  . ITH  INITIAL  DFL  AND  CURRENT 

BASE  PCINT 

00000620 

00000630 

c 

CALL  MERIT (PHI • SPI ) 

000006 A 0 
00000650 

IF(IPT.GE.O)  WRITF(6,707) 

IF ( I ThICE. EQ. 0 ) KE TURN 

00000660 

00000670 

l)ELsOEL  S 

C U T = 0 • 

00000680 

00000690 

c 

GO  TO  90 

00000700 

00000710 

c 

c 

SET  NEW  BASE  PCINT 

MAKE  PATTERN  MOVE 

00000720 

00000730 

c 

c 

EXPLORE  AKCJNi)  PATTERN 

00000740 

00000750 

20  3 

SSI  = S 

SSITST=F(SSI ) 

00000760 

00000770 

I T^R=  I TER  4-  I 

NPATM  = NPATV  4-  l 

. 

00000780 

00000790 

IF ( I TLR.GT.ITLIHGn  TC  700 

I F ( IPT.LT. 0 ) GO  TO  203 

00000800 

00000810 

WRITE (6,599)  I T F P 

WRITE( 6,595)  NPATM 

00000620 

00000830 

W R I T E ( 6 , 6 0 0 ) (PHI!  I ) , 1=  1 , N ) 

KRITE(6,60I)  SSI , DEL 

00000840 

00000850 

c 

c 

MAKE  PATTERN  MOVE 

00000860 

00000870 

c 

20  3 

DD  201  1=1, N 

00000880 

00000890 

THT ( I)=PSI  (I  ) 

PSI  ( I ) = PH I (I ) 

00000900 

00000910 

201 

PH(m=PHKI)  4 ALFAMPHKl)  - THT  ( I ) ) 
CALL  MERIT (PHI , SPI ) 

00000920 

00000930 

S=SPI 

I F ( IPT.NE . 1 ) GO  TO  202 

00000940 

00000950 

WRITE ( 6,606)  ( PHI  ( I ) , 1 = 1, N) 

W R I T F ( 6 , 60  1 ) SPI , DEL 

00000960 

00000970 

c 

202 

I C AL  L = 2 

00000980 

00000990 

c 

c 

MAKF  EXPLORATORY  MOVE 

00001000 

00001010 

c 

GO  TO  ISO 

00001020 

00001030 

c 

c 

IS  thl  present  value:  less  than  the  base 

POINT  VALUE 

00001040 

00001050 

250 

I F ( S.LT. S$  ITST ) GO  TO  200 

GO  T 9 100 

00001060 

00001070 

c 

c 

INTFRIJAL  SUBROUTINE  TC  MAKE  EXPLORATIONS 

ABOUT  PHI 

00001080 

00001090 

c 

150 

DO  180  K=l , N 

00001100 

00001110 

PH  I U LD=  PH  1 (K) 

STFPK0.05**>HIUL3 

00001120 

00001130 

1 F (ST EPK .E  0.0.0)  STEPK=0.05 

STfPK=SIGM (STEPK* DFL, XF LOCK) ) 

00001140 

00001150 

PHI  (K)  = PHICLD  4-  STEPK 

CALL  MCkI T (PHI , SP I ) 

00001160 

00001170 

lFMPT.ro.  1 UR  IT  F(  6,602)  ICALL,K,bPl,(PHl(l)  ,L  = 1,N) 
IF(S’M.LT.S)  GO  Tl  179 

00001180 

00001190 

XFI  » ( K ) = — XfLG(K) 

00001200 

56 


CALL  MERIT (PHI fSPI ) 
IF(LPT.ED.l)  WRlTt(6.602j 

) ICALL 

. * K • SPI  . 

(PHI ( L 1 .L= 

l.N) 

00001220 

OQOOlpTO 

IF(SPI.LT.S)  GfJ  TO  179 
PHIU)=PHICLD 

00001240 

00001250 

GO  TD  180 

179  S«SP1 

00001260 

00001270 

160  CONTINUE 

GC  TO  ( 160.260) . ICALL 

00001790 
0Q0Q1290 

000  013  00 


703  PSUII«PHI  (I) 

PT .GF . C) WRITE 


RETURN 

599  FORMAT ( * ♦»**  *.15) 


600  FORMAT ( • BASE  PT  % 1PM 7E1 5* 6#E 15.6, /, T9) ) 
606  FORMAT  ( 1 PATTERN  *.  1P4(7E15«6.E15.6./#T9)  ) 


601  FQKVAU6X,  •OBJ9,  2 Pc 2 5 . 6 , 5X r 1 DEL  *,E1  5.6) 

602  FORMAT!  IX,  212,  • QPJ  MPE14.6,'  TRIAL  ».4(6E14.6.  F14.6./.135 


609  FORMAT ( • 00 EL  •♦1PE15.6,*,  DELMIN  • t E 1 5.6t IX, // f * ITLIM  *fI6,*t 
1 I PT  1 1 1 3) 


607  FORMAT!9  TCTAL  NUMBER,  OF  NEW  BASE  POINTS  (ITERATIONS)  *,I5I 
701  FORMAT ( 9 OS  PARCH  T i rt M I NAT ED  BEC AUS F NUMBER  OF  ITERATION 


1LIMIT. • ) 

* 0 S EAP-C  H TERMINATED  BECAUSE  STEP  SIZE  LESS  THAN  LIMIT.*  ) 


707  FORM  AT ( • 0 STARTING  OVER  WITH  CURRENT  BASE  POINT  AND  INITIAL  DEL. * 
END  . 


00001340 

1 


00001360 


[WlWIiliU 


00001380 
" 01390 


00001400 
00014 


00001470 

00001430 


00001440 


00001460 


[ilium  mil 


) 000  0 1430 
00 


FUNCTION  A1(V) 


CnrtMGN/KLK/*/MD(K)  , S')(8)  f A I (8)*AJI25)»AK(B)  »U(50)  t 00000020 

L XMAXt  r.O 


00000020 


L4kl  UKkDl  fllUiUI  iW  WJ  VVliVI 


1 CLEM2C?Lt 9V15) 9ELEM3(2lv8915) vELEM4(?lv0f !5)9X(16!« I?  00000060 

!LU=I2+2  


A 2 =0 • 00000070 

W*  ( ( V-(  1)1  A >)/SiG)**2 


Y=RA(  lLU)-3.  00000090 

lcU.LE«5Q.)A2  = EXP(  V)»»Y  » EXPt-Q«5»W) 


RETURN  00000110 

CNO  


[flWIIlMWII 


2 C(  50)tkA(  I7),Ni,AN(8)  ,DI  A,  SIGtSl)LN(  81  fWMDLN(R)fTQ  00  300040 

CQMMuN  /MAN/  E1(S,15)»E2( 


00 


00000070 
0 


00000090 
0 


00000110 
00 


FUNCTION  A 3 ( V ) 

COMNION/MAOEP/I  l » J J t AL  F A l 8 f 5 0 ) , R ( 5 0 t 5 ) , N ( 50  ) t KC  ( 5 0 ) , A IF  AT  ( 50  ) 00000020 

1 tRK(*,5)  ,NC(50) »K,ERR{50)  ,EHNI?0) 000000 


COMMON /bLK /WMO ( 8 ) , SD(8)  , AI{8)fAJ(25)fAK(8),U<50)t  000 000 AO 

1 XMAX(50) >XMI  4(50) , RHGT (50) ,PHOX ,RHCH» API « SFFJ, 00000050 


2 C(50) #kA(17) 9Nl«AN(8)f0IAvSIG»SDLN(8) ,WMDLNI8)fTQ  00000060 

00000070 


X*(EXP<  V)  / AN  (3)  ) **2  000000 80 

ft  2=AN  ( l ) ♦ AN<  2)  / SORT  ( ( 1 ) **2  + A , » l AN  <4  ) »»2  ) *X  ) 30000000 


= R A ( 5 ) -3  • 00000100 

= B?*TQ  000001  Z 0 


TR  = ( < V- C r A > / S I G ) **2  00000120 

A3  = TS*EXP(-0.5*TRI  00000130 


RETURN  00000140 

END  00000130 


MILLER  DATA  SET  SO- I I I-  1 . -2 5 

C 82 10310 00 •500*  2000, 

.3138  .1368 

9 


.3198  .2421 

•4211  .1368 


. 3198.1368.3158 
. 3158.2421 


.3158  .2421 

.4211  .136 3 


.3158  .1368.4211 

•315°  .*421.3158 


,3 

.491 

.399 

.265 

.523 

— t. 

.3 

.481 

.323 

.228 

.444 

” * 

(ZERO  ADDITIVE) 

.4211  • 9 1 6 1 « 
.3158  .6891, 


.3158  .7971, 
.3158  .928  , 


1053  .6211, 

3158  .6921, 


3158  .7711, 

3153  .841  , 


.6171. 
613  . 


.6272.214 

.8812.291 


.6321.909 

.4711.704 


.7371.744 

.6801.774 


.6311.822 

.5001.604 


.6761.589 

6371.490 


.4211  .1368.3158 

3158.2421 


) .706 

.4491.169 

- 

L .561 

.407  .761 

- 

.3153 

.4211 


.315P.3158 

3158 


.5579 

.4526 


. 1055.1368 
4211.1368 


. 3158.1 56 


L 36 8 

L 368 

MILLER  DATA  SET  SD-IV-i 


C82 003 1000. 500.  2000. 

-3 


6231.050 


696  .916 
8 4 R .468 


424  .298 
6131 .037 


. 7041.671 


.5571.461 
, 3291.065 


.220  .397 
6551.533 


► 75 

-c 

) 89 

6021.250 
7 C 8 .568 
*6  12  .2  50 


,8081.861 
.3791.012 
,20*5  • 4 7H 


0P19O31OOO.  500.7000. 

*1 


3 . 1 


.3  .674  .812  .5221.327 


3 .462  .355  .261  .495 

.59?  -895  .5911-344 


.610  .869  .5651.310 


.502  .357  .263  . 


.450  .481  .352  .656 


.616  .462  .815 


MILLER  DATA  SET  S i)-V  I- 1 ..-2  5 (INCREASED  SOLIDS.  217,  24  MICRON  AL) 
0 5060310 00.  50 0.  2000 


.2957  .0986.2957  .799  .726  .4381.327 

3943 .2957 .800  .382 


2957  .3943  .645  .589  .390  .954 


.2957.29*7.0986  .457  .314  .229  .432 

943  .2957  .544  . 


MILLER  DATA  SET  SD-VII-1,-25  (187  24  MICRCN  AL  ♦ 17  FEO) 


i 

m 


iti*v  iivi  uviifwuQmnni 


.2957 
95 


.2957.0986 
957 


.29*7 

.9986 

. 

. 394  3 

9i 

.2957 


. 3943 
0986.2957 


.5391.303  .8421.778 

61 


mm  rwig.1 


,5571.146  .7401.601 
5311.039  .6871.435 


.762  .4951.132 


.3943 


.499  .615  .462  .921 
I . 


.515  .719  .479  .978 
5 


. 2957.0986 


.488  .762  .5331.048 
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